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We investigate the physical state of Hi absorbing gas at low redshift (z = 0.25) using 
a subset of cosmological, hydrodynamic simulations from the OWLS project, focusing in 
particular on broad (b m > 40 km s _1 ) H i Lya absorbers (BLAs), which are believed to originate 
in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which 
+3 includes radiative cooling by heavy elements and feedback by supernovae and active galactic 

nuclei, predicts that by z = 0.25 nearly 60 per cent of the gas mass ends up at densities and 
i ^ i temperatures characteristic of the WHIM and we find that half of this fraction is due to outflows. 

The standard H i observables (distribution of H i column densities N HI , distribution of Doppler 
Cn| parameters b Hl , b m — Nm correlation) and the BLA line number density predicted by our 

^ simulations are in remarkably good agreement with observations. 

BLAs arise in gas that is hotter, more highly ionised and more enriched than the gas 
giving rise to typical Lya forest absorbers. The majority of the BLAs arise in warm-hot 
(log (77 K) ~ 5) gas at low (log A < 1.5) overdensities. On average, thermal broadening ac- 
counts for at least 60 per cent of the BLA line width, which in turn can be used as a rough 
indicator of the thermal state of the gas. Detectable BLAs account for only a small fraction 
of the true baryon content of the WHIM at low redshift. In order to detect the bulk of the 
mass in this gas phase, a sensitivity at least one order of magnitude better than achieved by 
current ultraviolet spectrographs is required. We argue that BLAs mostly trace gas that has 
been shock-heated and enriched by outflows and that they therefore provide an important 
V window on a poorly understood feedback process. 
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1 INTRODUCTION 

The analysis of intervening H i Lyo- absorption in the spectra of dis- 
tant quasars (QSO) has become an extremely powerful tool to study 
the spatial distribution of the diffuse intergalactic medium (IGM) 
that follows the large-scale distribution of cosmological filaments, 
and to constrain the baryon content of the IGM as a function of 
redshift. At redshifts z > 3, more than 95 per cent of the baryonic 
matter resides in the form of photo-ionised, diffuse gas giving rise 
to the "Lyo- forest" in the spectra of distant QSOs (e.g. Rauch et al. 
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|1997[ l. As a consequence of expansion, the Lyo- forest thins out, and 
at z ~ the contribution of the Lyo- forest to the total baryon budget 
has decreased to ~ 20 per cent (e.g. |Penton et al.|2 004 Lehn er et al.| 
2007 ). At the same time, the formation of galactic structures and the 
gravitational heating of the IGM by collapsing large-scale filaments 
lead to a gradually increasing amount of shock-heated intergalactic 
gas at temperatures T > 10 5 K, which is referred to as the warm-hot 
intergalactic medium ( WHIM,|Cen & Ostriker|1999||Theuns et al] 

[T998l[Dave"et al.|2001||Bertone et al |2008| >. 

Since collisional ionisation determines the ionisation state of 
the shock-heated IGM, the neutral gas fraction in the WHIM is 
significantly lower, by at least one order of magnitude, than in the 
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photo-ionised IGM of the same density (e.g. Richter et al. 2008). 
Because of this very small neutral hydrogen fraction in the WHIM, 
most of the recent observational campaigns to study warm-hot in- 
tergalactic gas at low redshift have concentrated on intervening 
absorption by highly ionised metals in ultraviolet (UV) spectra of 
bright QSOs. In particular, five-times ionised oxygen (Ovi) has 
been used extensively to trace shock-heated intergalactic gas at low 
redshift and to constrain the baryon content of the WHIM (e,g. |Tripp| 
et al.|2000||Richter et al.|2004|[Dan"forth et al.|20Tffi|Danforth"&| 
Shull|2008; Thorn & Chen 2008b, ITripp et al. 2008 , Danforth et all 
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2010} . However, because O vi predominantly traces metal-enriched 
gas in a critical (in terms of ionisation balance) temperature regime 
at T as 3 x 10 K, and because the metals may well be poorly mixed 
on small scales I Schaye et al. 2007| the interpretation of intervening 
O vi absorbers is still contro versial (e.g. | Oppenheimer & Dave|2009l 



Tepper-Gar cia et aL]|2011[ |Smifh et al.|2011 1. In particular, it is 
not yet clear whether O vi absorbers predominantly arise in photo 



ionised (e.g. Thorn & Chen|2008a| > or collisionally ionised gas (e.g. 
Danforth & Shull 2008 ), or in complex absorbing structures with 
cool gas intermingled with warm-hot gas ( |Tripp et aT[ 2008 1. 

An alternative to highly ionised metals as tracers of warm-hot 
gas is offered by H i absorption. Due to the low neutral hydrogen frac- 
tion expected from collisional ionisation at temperatures T > 10 s K, 
Lyo- absorption from shock-heated WHIM filaments is expected to 
be very weak. In addition, H i absorption lines arising in gas at tem- 
peratures T > 10 5 K are expected to be relatively broad because of 
the effect of thermal broadening. Such broad (£> H i > 40kms~') and 
shallow ([Nhi/^hi] ~ 10"cnr 2 knr's or t (Hi) ~ 0.1) Lya absorp- 



tion features, the so-called Broad Lya Absorbers (BLAs; Richter 
|et al.|2006a[ l, are hence difficult to identify in the UV spectra of 
QSOs because of the limited signal-to-noise (S/N) and the low 
resolution of spectral data obtained with current space-based UV 
spectrographs. 

In spite of being observationally challenging, directly detecting 
the small amounts of neutral hydrogen in the WHIM in absorption is 
a feasible task. The first systematic studies of BLAs at low redshift 
have been conducted using high-resolution Hubble Space Telescope 
(HST) Space Telescope Imaging Spectrograph (STIS) spectra of 
bright QSO s <|Richter et al.|[2004{ |Sembach et al.|[2004] |Richter| 

Danforth et al. 



etal.|2006a, Williger et al. 2006; Lehn er et al.|2007| 



2010} . These studies indicate that BLAs may indeed account for 
a substantial fraction of the baryons in the WHIM at z ~ 0. They 
also show, however, that identification and interpretation of broad 
spectral features in UV spectra with limited data quality is afflic- 
ted with large systematic uncertainties. In particular, the effects of 
non-thermal broadening and unresolved velocity-structure in the 
lines lead to the occurrence of broad spectral features that do not 
necessarily arise in gas at high temperatures. The Cosmic Origins 
Spectrograph (COS; Green et a l.|2012} , a new UV spectrograph 
which has recently been installed on HST, is expected to substan- 
tially increase the number of BLA candidates at low redshift. Due 
to the limited spectral resolution of COS (~ 17kms~'), the system- 
atic uncertainties in identifying thermally broadened H i lines in the 
WHIM temperature range will nevertheless remain. 

To investigate the physical properties and spectral signatures 
of BLAs at low redshift, [Richter et aL] l |2006b) have studied broad 
H i absorption features using a cosmological simulation based on 
a grid-based adaptive mesh refinement (AMR) method ( |Norman| 
& Bryan 1999). Their simulation reproduces the observed BLA 
number density and supports the idea that BLAs trace (at least in 
a statistical sense) a substantial fraction of shock-heated gas in the 
WHIM at temperatures T ~ 10 5 - 10 6 K. However, since this (early) 



simulation ignored several important physical processes that are 
expected to affect the thermal state of this gas phase (i.e. energetic 
feedback, radiative heating and cooling by hydrogen and metals), 
it is important to re-assess the frequency and physical properties 
of BLAs using state-of-the-art cosmological simulations with more 
realistic gas physics. 

In this paper, we present a systematic study of BLAs at low 
redshift based on a set of cosmological simulations from the Over- 
whelmingly Large Simulations (OWLS) project (Schaye et al .|2010} . 
This work complements our previous study on intervening O vi ab- 
sorbers and their relation to the WHIM based on a slightly different 
set of OWLS simulations (Tep per-Garcia et al.|201l[ henceforth 
Paper I). The main features of the simulations we use are briefly 
described in Sec. [2] As we have done in Paper I for the case of low 
redshift O vi absorbers, we compare the predictions from our fidu- 
cial model to a set of standard H i observables, and discuss various 
physical properties of the general H i absorber population in Sec. [3] 
Given the dependence of the WHIM mass fraction predicted by sim- 
ulations on the particular implementation of the relevant physical 
processes reported in the past (e.g. Cen & Ostriker 2006), we invest- 
igate the impact of different physical models on the thermal state of 
the various gas phases in our simulations in Sec. [4] In this section 
we also present and discuss the results on the physical properties of 
the absorbing gas traced by BLAs. Finally, we summarise our main 
findings in Sec. [5] In the Appendix we include: a full description 
of our fitting algorithm (Appendix [A}; a detailed calculation of the 
observability of H i absorbing gas in terms of optical depth as a 
function of density and temperature (Appendix[B]l; a discussion of 
the convergence of our results with respect to the adopted physical 
model (Appendix |C), and with respect to the adopted resolution and 
simulation box size (Appendix [D^. 



2 SIMULATIONS 

The simulations used in this work are part of a large set of cos- 
mological simulations that together comprise the OWLS project, 
described in detail in Schaye et al. (2010, and references therein). 
Briefly, the simulations were performed with a significantly exten- 
ded version of the W-Body, Tree-PM, Smoothed Particle Hydro- 
dynamics (SPH) code gadget in - which is a modified version of 
gadget ii (last described in Springel 2005 1 -, a Lagrangian code 
used to calculate gravitational and hydrodynamic forces on a system 
of particles. The initial conditions were generated from an initial 
glass-like state ( White 1 996 ) with cmbfast (version 4.1; Seljak & Za- 



|ldarriaga|1996} and evolved to redshift z = 127 using the Zeldovich 
( 1970) approximation. 

The reference model, dubbed REF, in the OWLS 
framework adopts a flat ACDM cosmology character- 
ised by the set of parameters (O m , fit,. fi A , erg, n s , h] = 
{0.238,0.0418,0.762,0.74,0.95,0.731 as derived from the 
Wilkinson Microwave Anisotropy Probe (WMAP) 3-year dat^] 
( Sperg el et al.|2007|. This m odel includes star formation following 
|Schaye & Dalla Vecchia| ( [2008} , metal production and timed 
release of mass and heavy elements by intermediate mass stars, i.e. 
asymptotic giant-branch (AGB) stars and supernovae of Type la 
(SNIa), and by core-collapse supernovae (SNIIe) as described by 

1 These parameter values are largely consistent with the WMAP 7-year 
results iJarosik et al. 201 1 1, the largest difference being the value of <xg, 
which is 2 cr lower in the WMAP 3-year data than allowed by the WMAP 
7-year data. 
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Table 1. Overview of the simulations used in this study. All model variations 
are relative to REF. 



Model 


Description 


NOSN.NOZCOOL 


neglects SNe energy feedback and 




cooling assumes primordial abundances 


NOZCOOL 


cooling assumes primordial abundances 


REF 


OWLS reference model (see text for details) 


AGN 


includes feedback by AGN (fiducial model) 



(and other) models that are part of the OWLS project, see |Schaye| 
letaLlpOlO) , 



3 THE GENERAL H i ABSORBER POPULATION 

In this section we test the predictions of our fiducial model (AGN) 
against observations using a set of well-measured Hi observables: 
the H i column density distribution function (CDDF), the distribution 
of Hi line widths, and the correlation between Hi column density 
and line width. 



| Wiersma et al.|p009b| >. It further incorporates kinetic feedback by 
SNIIe based on the method of|Dalla Vecchia & Schaye] p508) , as 
well as thermal feedback by SNIa (Wiersma et al. 2009bJ. Radiative 
cooling by hydrogen, helium and heavy elements is included 
following the method of Wiers ma et al.| ( [2009a| l. The ionisation 
balance for each SPH particle is computed as a function of redshift, 
density, and temperature using pre-computed tables obtained with 
the photoionisation package cloudy (version 07.02.00 of the code 
last described by jFerland et al.||1998 1, assuming the gas to be 
optically thin and exposed to the |Haardt & M adau ( 2001 1 model for 
the X-Ray/UV background radiation from galaxies and quasars. It 
is worth noting that a simulation run that adopts the REF model, 
although with a slightly different set of values for the cosmological 
parameters (from WMAP7), has been shown to reproduce the Hi 
absorption observed at z = 3 in great detail ( Altay et al .|2011| >. 

Along with REF, we consider three further models from 
the OWLS suite respectively referred to as NOSNJVOZCOOL, 
NOZCOOL, and AGN. All these models differ from the reference 
model in one or more respects. NOSN _NOZCOOL neglects kinetic 
feedback by SNIIe, and the calculation of radiative cooling assumes 
primordial abundances. It is the most simple model in terms of input 
physics, and it is similar (and hence useful for comparison) to the 
simulation used by |Richter et"aL"l ( |2006b| . The model NOZCOOL 
assumes primordial abundances when computing radiative cool- 
ing, and the model AGA' includes feedback by active galactic nuclei 
(AGN) based on the model of black hole growth developed by |Booth| 
|& Schaye| ( [2009l see also |Springel et al.|2005) . 

All these simulations were run in a cubic box of 100/i~' co- 
moving Mpc on a side, containing 512 3 dark matter (DM) particles 
and equally many baryonic particles. The initial mass resolution is 
4.1 x 10* /r 1 M (DM) and 8.7 x 10V M Q (baryonic). The gravit- 
ational softening is set to 8 h' 1 co-moving kpc and is fixed at 2/i~' 
proper kpc below z = 3. 

In this study, we choose AGN as our fiducial model since it 
is the most complete model in terms of input physics. In addi- 
tion to reproducing various standard Hi statistics (see Appendix 
|C] >, this model has been shown to reproduce: the observed mass 
density in black holes at z = 0; the black hole scaling relations 
l |Booth & Schaye|2009) and their evolution ( |Booth & Schaye|2011) ; 
the observed optical and X-ray properties, stellar-mass fractions, 
star-formation rates (SFRs), stellar-age distributions and the ther- 
modynamic profiles of groups of galaxies (McC arthy et al.|2010| ); 
and the steep decline in the cosmic star formation rate below z = 2 
l |Schaye et al.|2010|[van de Voort et al.|201 Note that, while the 
Hi statistics predicted by the AGA' model are very similar to the 
predictions of the other models considered here (see Appendix [C^, 
there are notable differences in the temperatures of the gas traced 



by BLAs (see Fig 
in Sec.|44l Table fl 



Tip . We will address this point in more detail 
briefly summarises the relevant features of the 



models described above. For a more detailed description of these 



3.1 Synthetic spectra 

For a meaningful comparison to existing data, we generate 5000 ran- 
dom sightlines (1000 at five redshifts spanning the range < z < 0.5 
with step dz = 0.125) through the simulation box covering a total 
redshift path Az = 189, corresponding to an absorption path length 
A* = 275. 

We use the package specwizard written by Schaye, Booth, & 
Theuns to generate a synthetic spectrum for each sightline contain- 
ing absorption by H i Lyo- only. Briefly, we draw a random phys- 
ical sightline across the simulation box of size L, which is simply 
defined as the line between a given point on opposite faces of the 
box, and the collection of SPH particles with projected distances to 
this line smaller than their smoothing length. Next, we calculate the 
ionisation balance for each SPH particle as a function of redshift, 
density, and temperature, which we do using precomputed tables 
obtained with the photoionisation package cloudy (version 07.02 
of the code last described by Ferlan cTet al.|1998|, assumi ng the 
gas to be optically thin and exposed to the Haard t & Madau| ( |2001) 
model for the X-Ray/UV background radiation from galaxies and 
quasars. We divide the physical sightline into A p i x = [a(z) L/h] /Ax 
pixels of constant width Ax, where h and a(z) are the Hubble con- 
stant in units of lOOkms^'Mpc -1 and the expansion factor at the 
box's redshift z, respectively, and compute the smoothed ion density 
n ion , the ion density-weighted gas temperature, and the ion density- 
weighted peculiar velocity at each pixel. Proper distance bins of 
width Ax along the sightline are transformed into velocity bins of 
width Av = H(z)Ax, where H(z) is the Hubble parameter at redshift 
z; ion number densities are transformed into ion column densities 
via N ion = n ioll A.r, and gas temperatures into Doppler parameters 
using b T = ■\J2kT/m iol[ ., where k is Boltzmann's constant and m ioll 
is the ion's mass. The Hi optical depth t(u) at each pixel is com- 
puted assuming a thermal (i.e. Gaussian) profile, taking peculiar 
velocities into account, as described by jTheuns et al.|(l998| their 
Appendix 4). Finally, the optical depth spectrum is transformed into 
a continuum-normalised flux via F{v) — exp[-T(u)]. 

We convolve our spectra with a Gaussian line-spread function 
(LSF) with a full width at half-maximum (FWHM) = 7 km s" 1 and 
re-sample our spectra onto 3.5 km s -1 pixels. We add Gaussian noise 
to each spectrum assuming a flux dependent root-mean-square (rms) 
amplitude given by (S/NT 1 F(v), where S/N is the adopted signal-to- 
noise ratio. We assume a minimum, i.e. flux-independent noise level 
o~mm = 10~ 4 . This implies that our algorithm will underestimate the 
true column density of absorption features with a flux of the order 
of (or lower than) cr mm , which corresponds to a logarithmic central 
optical depth log r ~ 1 (see Appendix [A). Our choice of a perhaps 
unrealistically low value for cr m j n thus allows us to reduce the gap 
between the true and the fitted column density of saturated lines. 

We generate three sets of spectra, adopting S/N=10, S/N=30, 
and S/N=50, respectively. The spectra with S/N=10 and S/N=30 
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Figure 1. Column-density distribution function, /(Nm). obtained from 
observations and from 5000 simulated spectra with different S/N values 
spanning the redshift range [0, 0.5]. Error bars along the y-axis show Pois- 
son uncertainties computed using the tables by Gehrels ( 1986 1. Assuming 
/(Nm) k Nhi - '' for logarithmic column densities in the range [13.0, 15.2] 
([12.5, 15.2]), we find /? = 1.916 ± 0.044 (J3 = 1.917 ± 0.032) for S/N=10 
(S/N=30). The inset shows the comparison between the true (A-axis) and the 
fitted (y-axis) Hi column density integrated along each individual sightline. 
For clarity, only the result for S/N=10 is shown. The dashed line corresponds 
to a perfect match (see text for discussion). 



thus closely match the properties of the large sample thus far ob- 
tained with HST/STIS; these will be used in Sees. [3T2] [33] and[3g] 
to test the predicted Hi observables against observations; the syn- 
thetic spectra with S/N=50 are intended to investigate the physical 
properties of the H i absorbing gas following a statistical approach, 
in the remaining sections of the paper. 

Fitting of our 5000 synthetic spectra using the procedure de- 
scribed in Appendix|A|yields a total of 93430, 66705, and 28649 
components for S/N =50, 30, and 10, respectively. The resulting 
line-number densities and their corresponding Poisson uncertainties 
are (dJV/dz) = 494 ± 22 (S/N=50), 353 ± 19 (S/N=30), and 152 ± 12 
(S/N=10). For reference, the sample of 341 Lya absorbers at z S 0.4 
identified in seven FUSE+STIS spectra with average S/N > 10 by 
[Lehner et al. (2007) along an unblocked redshift path Az = 2.064 
yields (dN/dz) = 165 ± 13 at S/N « 10, which agrees (within the 
Poisson uncertainties) with our result at a similar S/N. 



3.2 Column-density distribution function 

In Fig.[T]we show the column-density distribution function (CDDF), 
/(N H i), obtained from our spectra with S/N=10 (red) and S/N=30 
(blue) spanning the redshift range [0, 0.5], together with results from 
different observations at similar redshifts using spectra with compar- 
able (average) S/N values. Assuming that the CDDF can be paramet- 
rized in the form of a single power-law, /(N HI ) oc N H i _/3 , we find 
P = 1.916 ± 0.044 (J3 = 1.917 ± 0.032) for S/N=10 (S/N=30) for 
logarithmic column densities in the range [13.0, 15.2] ([12.5, 15.2]). 
The lower limit in logN H i approximately corresponds in each case 
to the completeness limit as given by eq. (A3}, while the upper 
limit roughly indicates the column density above which our fitting 
algorithm underestimates the true Hi column density due to the 
minimum noise-level adopted (see Appendix |A}. 

The slope we obtain is in fairly good agreement with 
the slope measured from different observations. For z £ 0.4 
|Lehner et al.| ( |2007| their table 7) measure a range of values 



P = 1.52- 1.92 for absorbers in selected column-density inter- 
vals between log (Nni/cm -2 ) = 13.2 and 16.5, and line widths 
£>hi < 40 km s -1 or b a , < 80 km s _1 . If we extend the fitted column 
density range to log (N H i/ctrr 2 ) = 16.5, we find p = 1.90 ± 0.06 
and p = 1.95 ± 0.06 for S/N=10 and S/N=30, respectively. [Williger| 
|et al.| ( [20T0| l use a subsample from the Lehner et al. (2007} data and 
their own data at log (N H i/cnr 2 ) < 12.3, and find /? = 1.79 ± 0.1. 
|Dave et al] (2001) measure P = 2.04 ± 0.23 for absorbers with 
column densities log (N H i/ctrr 2 ) > 12.9 at a median redshift 
~Z = 0.17. Note, however, that a significantly shallower slope is found 
by |Penton et al.| l |2004} who identify 109 Lya absorbers at z < 0.069 
along 15 STIS spectra with S/N > 20, and measure p = 1.65 ± 0.07 
for logarithmic Hi column densities in the range [12.5, 14.5]. 

The amplitude of the CDDF resulting from the analysis of 
our synthetic spectra adopting different S/N is also in remarkable 
agreement with the observations. Note that the amplitude comes out 
naturally from our simulation, i.e. the CDDF has not been norm- 
alised to match the data in any way (even though that could have 
been justified because of uncertainties in the intensity of the UV 
background). At column-densities log (Nni/cm -2 ) < 15, our pre- 
dicted amplitude agrees well with the data all the way down to 
the lowest column densities measured, log (N Hl /cirr 2 ) = 12.3. At 
log (Nhi/cuT 2 ) > 15, the amplitude of our predicted CDDF appears 
slightly lower (or its slope is steeper) than the result by [Lehner] 
et al. ( 2007 1. Note, however, that their data point at highest meas- 
ured column-density bin has a rather large uncertainty. On the 
other hand, it is very likely that our choice of fitting parameters 
leads us to underestimate the amplitude of the predicted CDDF at 
log (NHi/cnr 2 ) > 14.5 by underestimating the true column dens- 
ity of saturated lines, as explained in Appendix [A] A comparison 
between the true and the fitted H i column densities integrated along 
each sightline reveals that our fitting procedure indeed yields integ- 
rated H i column densities which are systematically lower than the 
true total column density, in particular for log (N H i/ctrr 2 ) > 15 (see 
inset in Fig.^. This could explain the difference between our pre- 
dicted CDDF and the result by |Lehner et al.| l |2007^ at the high-N H , 
end. 



3.3 Line width distribution 

Fig.|2]shows the distribution of Doppler parameters, bm, obtained 
from our synthetic spectra for S/N=10 and 30 spanning the redshift 
range [0,0.5], together with the line width distributions obtained 
from data with comparable S/N values and redshi fts by|Lehner| 
et al. (2007 , green data points) and Danforth & Shull ( 2008 orange 
data points). The median values of our predicted distributions are 
b^a 30.4 km s -1 , 29.8 km s" 1 , and 29.4kms-» for S/N = 10, 30, 
and 50 (not shown), res pectively. All of the se agree well with the 
median value foun djyy Heap et al. (2002} , fo Hl = 27 kms"', by 



Shull et al. 



1 20001, fen i = 28 km s 1 , and with the median value 
b m = 31 kms" 1 for the full[] 



Lehner et al. 



(2007 



sample. Note that 
all of these values are signi ficantly larger than th e median value 



&hi - 21 km s measured by Dave & Tripp 



1 2001 ). Our simulation 



shows a lower fraction of broad (Z? H l > 40 km s~ 1 ) absorbers when 
compared to the Lehner et al. (2007 1 b- value distribution, but our 
results compare well to the line width distribution from Dan forth "&| 
|Shim](2008} . 

The predicted median /? H i -values indicate that a lower S/N 
value systematically shifts the line width distribution to slightly 
larger values. Yet, the number of components with £> HI > 40kms~' 
relative to the total number of components identified in each case 
decreases from ~ 26 to ~ 23 per cent when the adopted S/N value 
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Figure 2. Distribution of bu r values from observations and from 5000 simu- 
lated spectra with S/N=10 (red), and S/N=30 (blue), spanning the redshift 
range [0, 0.5]. Data points from Lehner et al. 1 2007 green) and |Danforth| 
|& Shull|j2008| orange) with y-error bars showing Poisson uncertainties. 
The dashed vertical line indicates our adopted minimum allowed b- value. 
The distributions from both data and simulations have been binned using 
A&hi = 5 km s . 



decreases from 50 to 10. Here, two competing mechanisms are at 
work: On the one hand, a low S/N value results in a stronger blending 
of narrow components into (artificial) broad features. On the other 
hand, since broader lines are shallower (at a given column density), 
and thus more difficult to detect at low S/N, the number of broad 
components detected decreases with decreasing S/N. Compared to 
a higher S/N value, the net effect of a low S/N value is to yield a 
smaller number (both relative and absolute) of broad absorption 
features (at a given resolution and sensitivity). 



3.4 The &h i - N H i distribution 

Last, we compare the £>hi — Nhi distribution obtained from our 
simulated spectra with S/N=30 and S/N=10 to two different sets 
of observations used for the comparison of our predicted CDDF 
and the line width distribution discussed in the last sections. To 
this end, we bin the lines from observations and from our synthetic 
spectra in N HI using A log(N H i/ctrr 2 ) = 0.3, and compute the 
median £>Hi-value, and 25-/75-percentiles in each bin. The result 
is shown in Fig. [3] The & H i _ Nhi distribution from our simulated 
spectra matches the observations well within the uncertainties. 
Even the drop in &hi observed at low Nhi in the |Dave & Tripp] 
< |200ip is well reproduced by our simulation. Note that lines with 
log(N H i/cnT 2 ) < 13.4 identified in spectra with S/N=30 generally 
have larger widths. This is a consequence of the fact that, at a fixed 
column density, lines with a given width are shallower with respect 
to narrower lines, and they can only be detected if the S/N is high 
enough. 

Summarising, we conclude that the H i observables predicted 
by our fiducial model are in excellent agreement with observations. 
This agreement may be surprising in view of the uncertainty in the 
input physics used in our simulation. However, in Appendix |C| we 
show that these results are quite robust against the model variations 
with respect to our fiducial model considered here (see Sec. |4.1[ l. We 
now proceed with the analysis of the physical conditions in low-z 
H i absorbers. 




Figure 3. Z>hi - Nhi distribution obtained from 5000 simulated spectra with 
S/N=10 (red), and S/N=30 (blue), spanning the redshift range [0, 0.5], and 
measurements by Lehner et al. 1 2007 green data points from their Table 3 for 
^Hi > km s~') and jDave & Trippj |200l| olive data points). Points show the 
median b\\ r value in each bin of size A log(Nn i / cirT 2 ) = 0.3, while the error 
bars parallel the y-axis correspond to 25 and 75 percentiles, respectively. 
Note that the red (blue) histogram on the right sub-panel corresponds to the 
red (blue) histogram in Fig. [2] but with a different binning and on a different 
scale. 

3.5 Physical state of the H i absorbing gas 

In this section we present and discuss the physical properties of 
the gas detected via H i absorption in our fiducial model (AGN; see 
Tab.[T|. The method we use is similar to the method described in 
Paper I, in which we used optical-depth weighted quantities. Briefly, 
to compute the desired Hi optical-depth weighted quantity (e.g. 
density) associated with a given absorption line, we first compute the 
optical-depth weighted density in redshift space along the sightline 
as in |Schaye et al.|fl999) . Next, we compute the average of the 
optical-depth weighted density over the line profile, weighted again 
by the optical depth in each pixel and assign this last weighted 
average to the line. In concordance with Paper I, in the following 
we shall denote quantities weighted by H i optical depth by adding 
a corresponding subscript; thus, for example, the H i optical-depth 
weighted temperature is denoted by T SI . We refer the reader to 
Sec. 5. 1 of Paper I for a more detailed description about our method 
for computing optical depth-weighted quantities. 

For simplicity, we obtain a new line sample of H i absorbers 
identified in synthetic spectra with S/N=50 generated from 5000 
sightlines across a simulation box at a single redshift^] z = 0.25, 
spanning a total redshift path Az = 187.5, which corresponds to an 
absorption path length Ax = 270. These spectra have been fitted 
following the method described in Appendix [A] 

We restrict our analysis to "simple", i.e. single-component, 
absorbers, unless stated otherwise, We define an absorber i as 
'simple' if the velocity distance from its centre to any other 
component j along the same sightline satisfies Av > 2<r b , where 
al = 0.5[fo Hl (() + £> HI 0)]. Absorption lines that do not satisfy this 
condition are referred to as 'complex'. 

2 Note that our chosen redshift is slightly higher than the median redshift of 
most Hi absorption-line studies at low redshift (e.g. ~z& 0.17 in |Lehner et al. | 
2007) . Although some evolution does take place from z = 0.5 — * 0.0, we do 
not expect the choice of this particular redshift to affect our conclusions in 
any significant way. 
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Table 2. Line-number density, (dN/dz), total baryon content in Hi, fim, 
and total baryon content in gas traced by H i, ilb(H i), related to simple H i 
absorbers identified in 5000 spectra with S/N=50, 30, and 10 at z = 0.25. 





S/N=50 


S/N=30 


S/N=10 


components (rel. to S/N=50) 


1 


0.72 


0.31 


simple absorbers (rel. to total) 


0.47 


0.53 


0.65 


(dN/dz) a 


460 ± 22 


332 ± 18 


144 ± 12 


a m [io- 7 ] b 


1.20 


1.15 


1.12 


n b (Hi)/£j b 


0.57 


0.47 


0.29 



a Quoted uncertainties are purely Poissonian. For comparison, Leh ner et al.| 
(2007) obtain (dN/dz) = 165 ± 13 at S/N a 10. 

b Total baryon content in H i obtained by adding the column densities of 
all identified H i components. The true total baryon content in H i along the 
fitted sightlines atz = 0.25 is £l H i = 2.11 X 10~ 5 . 

Table[2]contains various statistical and physical quantities res- 
ulting from the analysis of these new line sample, such as the relative 
number of identified components, the relative number of simple ab- 
sorbers, the line-number density, (dN/dz), the total baryon contenj^] 
in H i, fl H i. an d the total baryon content in gas traced by H i, Qt,(H i) 
(see also Sec. [43}. Note that the statistical and physical properties 
of this new sample are very similar to the corresponding properties 
of the line sample discussed in Secs.|3.2|-|3.4| 



3.5. 1 Physical density and absorber strength 



As previously noted by several studies (e.g. Schaye et al. 1999, Dave 
|et al.| [T999), there exists a tight correlation between Hi column- 
density, Nhi, and overdensit}]^] A = p h / (pb), of the absorbing gas 
usually parametrized in the form of a power-law, A/A = (Nm/No)". 
Due to variations in the (local) ionising radiation field, the influence 
of other heating mechanism (shocks), and other factors such as 
the geometry of the absorbing structures, etc., this relation has an 
intrinsic scatter, which decreases with increasing redshift ( |Dave et al"| 
[1999). 

The relation between overdensity and H i column density for 
the diffuse IGM has been derived analytically by Sc haye| ( |2001| l, 
who assuming local hydrostatic equilibriurr^jand optically thin gas 
finds 



log A oc 



3 + a (1 - 2b) 



log N H , 



log(l+z) 



(1) 



In the above equation, a is the slope of the temperature-density 
relation, T = T A", which results from the balance between photo- 
heating and adiabatic cooling ( |Hui & Gne din 19971, and —b is 
the power of the temperature in the expression for the Hi recom- 
bination rate coefficient which behaves as oc T~ b . If re-ionisation 
of the IGM takes place at sufficiently high redshifts, its imprints 
on the thermal state of the IGM are eventually washed out, and 
the slope of temperature-density relation is expected to reach 



The total baryon content in H i is computed via 



Ohi 



4 The 

(Pb) = 



model 



baryonic density 
(/1/0.73) 2 (1 + z) 3 g cm- 
5 The assumption of 'local hydrostatic equilibrium' implies that the size of 
a self-gravitating gas cloud is of the order of the local Jeans length. 
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Figure 4. Ah i - Nh i relation for simple H i absorbers identified in spectra 
with different S/N values at z = 0.25 for. Points show the median overdensity 
in each bin of size Alog(NHi/cnT 2 ) = 0.1. Error bars along the y-direction 
correspond to 25 and 75 percentiles in each bin. For clarity, only the error 
bars for the S/N=5() result are shown, but they are similar for S/N=30 and 
S/N=10. The horizontal arrows indicate the formal completeness limit for 
each adopted S/N value as given by eq. \A3) . Note the deviation of the 
Ah i - Nhi relation from a single power-law for column densities below 
the formal completeness limit (at a given S/N) and in the column density 
range for which the H i hya line generally saturates (shaded area). A power- 
law with the theoretical expected slope 0.738 (but with arbitrary amplitude; 
|Schaye|2 001 1 has been included to guide the eye (black, dashed line). 



an asymptotic limit determined by the temperature dependence 
of the Hi recombination rate. More specifically, at low redshift 
a -> 1/(1 + b). We fintj^ = 0.755 ± 0.001 in the temperature range 
[10 3 , 5 x 10 4 ] K, and hence a — > 0.570. Inserting this value into 
eq. |Q} gives A oc N H i°- 738 ■ (1 + z)~ 3 36 . Thus, the value of the amp- 
litude A in the A - N H r relation decreases with redshift, implying 
that absorbers of a given column density trace gas at higher over- 
densities at lower redshift. 

In Fig. [4] we show the Am - N H i relation resulting from our 
fiducial model for all simple Hi absorbers identified in spectra 
with different S/N values at z = 0.25. The turn-over in overdensity 
at column densities below the sensitivity limit (eq. |A3[ > for each 
adopted S/N is caused by errors in the measured N H i below this 
limit. Note also the deviation of the Ahi - Nhi relation from a single 
power-law at column densities log(N H i/cirr 2 ) > 14.5 (indicated by 
the shaded area), which corresponds approximately to the column 
density for which the Lya line saturates. This is in part due to the 
inability of our algorithm to properly fit saturated lines. 

Performing a least-square, error-weighted fit to the S/N=50 
result for column densities above the corre sponding sensitiv- 
ity limit (log(N H i/cm -2 ) = 12.3; see eq. 
log(N H ,/cirr 2 ) < 14.5, we find a = 0.786 : 



A3 i 



and restricted to 
: 0.010 and A = 48.3, 
normalised to No = 10 14 cirr 2 . The resulting slope for S/N=30 
(S/N=10) is a = 0.798 ± 0.010 (a = 0.847 ± 0.014) in the logar- 
ithmic column density range [12.5, 14.5] ([13.0, 14.5]), where the 
lower column density limit is given by eq. ( |A3[ >. A power-law with 
the theoretical expected slope 0.738 and arbitrary amplitude has 
been included in this figure for reference (black, dashed line). 

For comparison, Da ve et al"] ( |2010| their equation 3) find 



We compute the recombination rate coefficient for recombination case A 
numerically, and fit a power law in the given temperature range. 
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Figure 5. Ratio of H i line thermal width, (6 r / km s~ 1 ) = 1 2.9 V7hi/10 4 K, 
to total line width, bm, as a function of the optical depth at line centre, 
ro, inferred from Voigt profile fits, for simple Hi absorbers identified in 
spectra at z = 0.25 adopting different S/N values. The top axis indicates 
the corresponding Hi column density assuming £>hi = 30 tans -1 , which 
corresponds to the median value of the Doppler parameter distribution (see 
Sec. |3.3) . The points show the median (frr/^Hi) value in each bin of size 
A log ro = 0.25 indicated by the error bars parallel to the .v-axis; the lower 
and upper error bars parallel to the y-axis correspond to the 25 and 75 
percentiles in each bin, respectively. For clarity, only the error bars for the 
S/N=50 result are shown, but they are similar for S/N=30 and S/N=l(). 
Note that the values corresponding to S/N=30 (blue crosses) have been 
slightly shifted for display purposes. The horizontal arrows indicate the 
formal completeness limit for each adopted S/N value as given by eq. {A3}. 
The dashed horizontal lines enclose the range 0.5&hi < fer < bui, and have 
been included to guide the eye. The gray, shaded area indicates the central 
optical depth (or column density) range for which the H i Lya line generally 
saturates. The hatched area indicates the typical range in ro for BLAs, 
- 1 .34 < log tq < 0, detected in spectra with S/N=50 (see Sec. |4.2.1) . 

a = 0.741 ± 0.003 and A = 38.9 at z = 0.25, for absorbers arising 
in gas with temperatures log (77 K) < 4.5 in their simulation. 
If we restrict our sample to single-component absorbers with 
log(7WK) < 4.5, we find a = 0.786 ±0.014 and A = 45.1. We 
note that we do not rescale the amplitude of the UV background in 
our simulation, while |Dave et al.| ( |2010[ l adjust its amplitude by a 
factor 3/2 to bring their predicted evolution of the H i optical depth 
into better agreement with observations. 

3.5.2 Gas temperature and line width 

A matter of interest is to which extent the measured Hi line 
width can be used to estimate the temperature of the Hi absorb- 
ing gas. We explore this by comparing the Hi thermal line width, 
b T = 12.9 -^7hi/10 4 K, computed from the optical depth-weighted 
gas temperature, 7hi, to the total Hi line width, bm, as a function 
of the line strength as given by the optical depth at the line centre, 
To. The optical depth at the line centre is computed using the Hi 
column density and the H i line width inferred from a Voigt profile 
fit to the line (see Appendix|A]l. We bin the ratio (b T /&hi) in T o> an d 
plot in Fig.|5]the median value and the 25-/75-percentiles in each 
bin as a function of To, for all single component absorbers identified 
in our synthetic spectra with different S/N values at z = 0.25. 

We see that thermal broadening becomes increasingly import- 
ant with increasing line strength (or Hi column density), and it 
contributes with at least 50 per cent to the total line width, i.e. 
bj > 0.5 bni, irrespective of the line strength and the adopted S/N 



value. The temperature of gas giving rise to absorption lines with 
central optical depths in the range 1 < To < 8.5 (corresponding to 
strong lines) on average contributes with at least 90 per cent to the 
total line width, i.e. b T > 0.9 b H] . The lower {b T lb m ) value shown 
by highly saturated lines, i.e. lines with To > 8.5 (gray, shaded area) 
is most probably due to the uncertainty in the fit parameters of 
such lines. These results suggest, in view of the tight Am - N H i 
relation discussed above, that absorption arising in low density gas 
is subject to more significant non-thermal (i.e. Hubble) broadening 
than gas at higher density. This is consistent with the idea that low 
density, unbound gaseous structures are subject to the universal 
expansion, while gas at higher densities residing closer to galax- 
ies may have detached from the overall expansion. We find (not 
shown) that the line width correlates well with the gas temperature 
for log (N H i/cm _2 j > 13, but that it is a poor indicator of the thermal 
state of the gas for lower column densities. 

Fig. [3] shows also that lines with central optical depths cor- 
responding to H i column densities below the formal completeness 
limit for each adopted S/N value (indicated by the arrows) as given 
by eq. {A3} have, on average, b T > bm, which is un-physical. These 
lines correspond to absorption by gas at high temperatures, which 
gives rise to very shallow and extremely broad features that are (in- 
correctly) fitted with several components, thus yielding line widths 
that are narrower than allowed by the gas temperature. 

In the central optical depth range typical for BLAs detected in 
spectra with S/N=50, -1.34 < logT < (see Sec. |4.2.1| > indicated 
by the hatched area, the contribution of thermal broadening to the 
line width amounts to 60 to 90 per cent. If non-thermal processes 
(e.g. turbulence) contribute to the line broadening in such a way 
that the total line width is given by = b\ + b\ v where b nt is the 
non-thermal broadening (as would be the case for a purely Gaussian 
turbulence field), then the ratio of non-thermal broadening to total 
line width can be important, even though the thermal contribution 
is substantial. Take as an example the maximum, average thermal 
broadening to total line ratio for BLAs b T /bm = 0.9; this value 
together with £> m = b\ + b\ t implies b^/b^ = 0.4. 

3.5.3 The (N# ; , b H ,)-plane 

A deeper insight into the physical state of gas giving rise to Hi 
absorption identified in real QSO spectra can be gained by study- 
ing the relation between selected physical quantities and the line 
observables, N H[ and fc H i, simultaneously. We have followed such 
an approach in Paper I in order to study the physical conditions of 
O vi absorbing gas, and now apply it to study gas traced via H i ab- 
sorption. We focus on four quantities: gas temperature T Hu neutral 
hydrogen fraction (nm/n^m, total hydrogen column density Nh, 
and gas metallicity Z. Note that gas temperature and metallicity 
are 'true' optical depth-weighted quantities, while total hydrogen 
column-density and ionisation fraction are 'derived' quantities. For 
instance, the neutral fraction is computed using the optical depth- 
weighted hydrogen particle density, (n H )ai> and the optical depth- 
weighted gas temperature, using pre-computed tables obtained with 
the photoionisation package cloudy (version 07.02 of the code last 
described by Ferland et al. 19981, as described in Sec. [2] 

For each of these physical quantities we proceed as follows: 
First we compute the desired physical quantity, e.g. 7hi» for each of 
the simple absorbers in our S/N=50 line sample. We then divide the 
(Nhi, &Hi)-plane into cells, and compute the median value for the 
desired quantity using the values of all absorbers with (N H i, £m) 
values in that cell. Fig. [6] displays the result for temperature (top- 
left), neutral hydrogen fraction (top-right), metallicity (bottom-left), 
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Figure 6. Temperature (top-left), neutral hydrogen fraction (top-right), metallicity (bottom-left), and total hydrogen column density Nh (bottom-right) of the gas 
at z = 0.25 traced by simple Hi absorbers identified in spectra with S/N=50, as a function of the line parameters Nhi and bm. The colour coding shows the 
median value in each cell of the corresponding physical quantity. The white contours show the distribution (by number) of the sample of simple H i absorbers, and 
they are identical in all panels. These contours enclose, starting from the innermost, 20, 40, 60, 80, and 90 per cent of the total number of single-component Hi 
absorbers. The horizontal and diagonal dashed lines indicate the usual BLA selection criteria bm > 40kms~' and (Nhi/^hi) = 3 x 10 12 cirr 2 kirr's (S/N) -1 
(e.g. |Richter et al.|2006a) , respectively. Absorbers within this boundaries typically trace gas at log(77 K) > 4.7, have total hydrogen column densities 
18 < log(NH/cirT 2 ) <, 18.9, and low ionisation fractions log(nHi/nH)m i> -5.5. Also, the (local) metallicity of the gas traced by these absorbers is typically 
log(Z Hl /Z G ) > -2.5. 



and total hydrogen column density (bottom-right). The colour code 
indicates the median value of the corresponding physical quant- 
ity. For reference, we include contours (white solid curves) show- 
ing the distribution by number of the simple H i absorbers on the 
(Nhi, &m)-plane. These contours enclose, starting from the inner- 
most, 20, 40, 60, 80, and 90 per cent of the total number of H i com- 
ponents. The dashed horizontal and diagonal lines at the top-right 
corner of each panel indicate, respectively, the BLA selection cri- 
teria b m > 40kms~ 1 and (N H i/fc H i) = 3x lO^cm^knr's (S/N)" 1 
JRichter et al.|2006aj adopting S/N=50. 

There are several interesting features in this figure. First, all 
four physical quantities appear to have a relatively simple depend- 
ence on Nhi and bm- The temperature of the gas (top-left panel), 
for example, shows a positive correlation with the line width, which 
appears to be tighter for absorbers at a given Nhi (range), in agree- 
ment with the results presented in Sec. |3.5.2] In this respect, note 
the population of narrow (£> H i ~ lOkms 1 ), low-column density 
(N H i < 10 12 cnr 2 ) absorbers with high (log (7W K) > 4.8) optical 
depth- weighted temperatures. These correspond to the absorbers 



with column densities below the formal completeness limit and with 
br > bm, previously discussed. 

The neutral hydrogen fraction (top-right panel) increases with 
N H i, but strongly decreases with bm- This can be interpreted as 
a temperature-dependence, given the positive correlation between 
7hi and bm- Correspondingly, the total hydrogen column density 
(bottom-right panel) increases with both Nhi and bm- The bottom- 
left panel shows that the (local) metallicity of the gas is strongly cor- 
related with the line width. Given the correlation between gas tem- 
perature and line width shown the top-left panel, this suggests that 
there is a correlation between gas temperature and (local) gas metalli- 
city. This correlation is very likely a consequence of strong feedback. 
Indeed, high-temperature, high-metallicity absorbers could be tra- 
cing shock-heated, enriched outflows in the surroundings of galaxies 
that have not had enough time to mix with the surrounding gas and 
to cool down, whereas low-temperature, low-metallicity absorbers 
could be tracing both gas that has not yet been impacted by outflows, 
and wind material that has been ejected at redshifts high enough for 
it to cool down and to dilute its metal content in the ambient gas. 

The BLA selection regime defined by the dashed lines in 
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Table 3. Definition of the various gas phases we consider in terms of temper- 
ature and density thresholds. For example, we define the warm-hot, diffuse 
gas (WHIM) to have temperatures T ^ 5 X 10 4 K and overdensities A < 10 2 . 





Overdensity (A) 


Temperature (T[ K]) 


cool 




< 5 X 10 4 


warm-hot 




> 5 X 10 4 


diffuse 


< 10 2 




condensed 


> 10 2 




star-forming 


> 3 x 10 5 


(EoS) a 



a We consider 'star-forming' the gas with physical densities that exceed 
our adopted star-formation threshold njj = 0. 1 cnT 3 and which is allowed to 
form stars. The temperature of this gas phase is set by an imposed equation 
of state (EoS) of the form P oc p 4 ' 3 . This gas phase can be thought of as 
the inter-stellar medium (ISM). Note that, although this gas is cold, it is not 
included in the gas phase defined as 'cool' . 

each panel reveals a population of Hi absorbers tracing highly 

ionised (log(n H i/n H )Hi 6) gas with median temperatures 

log (7^7^) <: 4.7, median (local) metallicities \og(Zm/Z ) > -2.5, 
and total hydrogen column densities log(N H /cm -2 ) a 18.7, which is 
almost an order of magnitude higher than the total hydrogen column 
density of typical Lya forest absorbers (see also Fig.|12|l. According 
to our previous interpretation of the r H i -Z H i correlation, these 
results suggests that BLAs may be tracing galactic outflows. We 
will come back to this point in Sec. |4.4.2"] 

The results presented in this section indicate that the H i column 
density of unsaturated absorbers is a reliable tracer of the underlying 
physical density of the gas giving rise to the detected H i Lya absorp- 
tion. Moreover, the temperature of the absorbing gas may be roughly 
estimated from the measured line width, as suggested by the average 
contribution of thermal broadening to the total line width of these 
absorbers. Finally, H i absorbers subject to the commonly adopted 
BLA selection criteria trace gas which appears to be physically 
distinct from the gas traced by typical Lya forest absorbers. 



4 THE WARM-HOT DIFFUSE GAS 

In the next sections we explore in detail the effect of feedback and 
metal-line cooling on the physical state (i.e. density and temperature) 
of the gas in our simulations. Also, we investigate the H i absorption 
signatures of warm-hot diffuse gas, and the physical properties of 
the gas traced by broad H i absorption features (BLAs) identified 
in synthetic QSO absorption spectra. For this purpose, we use the 
H i sample from our fiducial model AGN presented in Sec. |3.5| and 
generate similar samples for all the other model runs. Comparison 
between model predictions and observations (whenever possible) 
are done exclusively for our fiducial run. 

4.1 Model-dependence of the predicted warm-hot gas mass 

On super-galactic scales, two mechanisms are able to shock-heat 
intergalactic gas to temperatures T > 5 x 10 4 K: a) galactic out- 
flows driven by SNII explosions and by AGN activity; b) accretion 
shocks caused by infall onto the potential wells of dark matter 
halos. We have selected four model runs from the OWLS project, 
NOSNJNOZCOOL, NOZCOOL, REF, and AGN, to investigate the 
effect of each of these mechanisms on the properties of diffuse gas 



and its imprints on simulated absorption spectra. Note that these 
models have already been described in Sec. [2] (see also Tab.[TJ. 

We are interested in the predicted distribution of gas mass 
among the various (gas) phases -in particular the warm-hot dif- 
fuse phase- defined in Tab. [3] We adopt a temperature threshold 
at log (77 K) = 4.7 (or T ~ 5 X 10 4 K) and a density threshold at 
A = 10 2 to distinguish these gas phases. The density threshold has 
been chosen so as to roughly separate unbound gas from collapsed 
structures (at z = 0.25). The temperature threshold is motivated by 
the bi-modality in the gas mass distribution predicted by the models 
considered here (see below). Note that our value is somewhat below 
the 'canonical' but to some extent arbitrary T = 10 5 K commonly 
adopted to distinguish cool from warm-hot intergalactic gas (but see 
|Wiersma et al.|2010| >. 

The distribution of gas mass as a function of temperature and 
(over-)density predicted by the different models is presented in Fig.|7] 
The coloured areas show the cumulative gas mass (in per cent) in- 
dicated by the colour bar on the right. The vertical (horizontal) 
solid line in each panel indicates the density (temperature) threshold 
separating the various phases. Star-forming gas, which is defined 
as gas with physical densities njj > 0.1 cirT 3 (oiQa a 3 X 10 5 at 
z = 0.25), is shown to the right of the blue, vertical dashed line in 
each panel. The percentage included in each separate region indic- 
ates the baryonic mass in the corresponding (gas) phase relative to 
the total baryonic mass. In particular, the number (orange) in the 
top-left corner of each panel gives the mass fraction of gas with 
r>5x 10 4 K and A < 10 2 , i.e. warm-hot diffuse gas. The starred 
percentage indicates in each case the baryonic mass confined in stars. 
Note that our adopted temperature threshold seems appropriate to 
separate cool, photo-ionised from shock-heated gas; the gas mass 
distribution at A < 10 2 is clearly bimodal, with two phases having 
significant gas mass fractions above and below log (77 K) = 4.7. 
The dotted contours indicate the neutral hydrogen fraction, (nni/nn), 
as function of density and temperature at z = 0.25; it has been com- 
puted using pre-computed cloudy tables as described in Sec. [2] The 
logarithmic (nHi/ n H)- value is indicated next to the corresponding 
contour in the top-left panel. 

The sequence of models given by moving clock-wise from 
the top-left panel is essentially a sequence of increasing feedback 
strength (and model complexity). The mass fraction in warm-hot 
diffuse gas in the model NOSNJJOZCOOL (top-left panel) indicates 
that by z = 0.25 roughly 30 per cent of the gas mass is shock-heated 
to temperatures log (77 K) > 4.7 by gravity alone. In the absence of 
any feedback on galactic scales, a large fraction of the gas that is 
accreted via gravitational infall at higher redshifts is able to cool 
and fuel star formation, with nearly 15 per cent of the gas mass 
ending up in stars by z = 0.25. The cool, photo-ionised diffuse gas 
at log (77 K) < 4.7 and A < 10 2 in this model contains roughly 40 
per cent of the total gas mass. 

Moving on to the top-right panel (model NOZCOOL), we see 
that nearly 45 per cent of the baryonic mass in the simulation is in 
the form of warm-hot diffuse gas as a consequence of the kinetic 
energy released by supernova explosions. This corresponds to an 
absolute increase in mass of 15 per cent in this gas phase compared 
to NOSNJJOZCOOL. At the same time, the mass fraction in the 
cool diffuse IGM predicted by NOZCOOL decreases with respect 
to NOSN_NOZCOOL, but by a far smaller amount (~ 4 per cent in 
absolute terms). Thus, about 10 per cent of the gas mass that ends 



7 The relation between hydrogen particle density nn and (baryonic) over- 
density A is given by eq. {57} in AppendixlB] 
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Figure 7. Distribution of gas mass over various phases predicted by different models, indicated on the top-right corner of each panel (see Tab.[TJ. The vertical 
(horizontal) dotted line indicates the density (temperature) threshold at A = 10 2 (T = 5 X 10 4 K) that separates unbound (cool) from collapsed (warm-hot) gas. 
The region to the right of the blue dashed line shows the high-density, star-forming gas with physical densities nj^ ^ 0.1 cirT 3 (or A > 3 X 10 5 at z = 0.25). The 
coloured areas show the cumulative gas mass (in per cent) indicated by the colour bar on the right. The percentages in each panel show the baryon mass fraction 
in the corresponding phase defined by the temperature/density thresholds. We have highlighted the baryon mass fraction in the diffuse, warm-hot gas in orange. 
The starred percentage indicates in each case the baryonic mass fraction in stars. The dotted contours, which are identical in all panels, give the neutral hydrogen 
fraction (nHi/nn), of gas in ionisation equilibrium as a function of A and T; the logarithmic (nHi/nn)-value is indicated next to the corresponding contour in the 
top-left panel. The left panels show the effect of our two most extreme scenarios, i.e. including both feedback by SNII and AGN with respect to neglecting 
feedback altogether; the right panels show the effect of neglecting radiative cooling by heavy elements. Clearly, feedback by SNII (and AGN) heats a significant 
amount of gas above temperatures T = 5 X 10 4 K, with the WHIM fraction increasing from 28.5 per cent (top-left) to 58.8 per cent (bottom-left). Interestingly, 
the IGM fraction only changes by ~ 10 per cent, indicating that feedback shifts a large fraction of the ISM from haloes into intergalactic space. 



up in the WHIM by z = 0.25 must be removed from a gas phase 
other than the cool diffuse IGM. The significantly lower mass in 
stars in the NOSNJNOZCOOL model compared to the NOZCOOL 
model strongly suggests that at higher redshifts supernova feedback 
shock-heats and blows a significant fraction of the ISM out of halos, 
which ends up in the WHIM by z = 0.25. 

Including radiative cooling by heavy elements (model REF; 
bottom-right panel) has a negligible effect on the WHIM and the 
cool diffuse IGM, suggesting that metal-line cooling in the these gas 
phases is inefficient, either because their density is low, or the metals 
they contain are not yet well mixed, or perhaps because the level 
of enrichment is low, or a combination of them all. Interestingly, 
the REF model predicts a much higher mass fraction in stars and 
star-forming gas with respect to NOZCOOL, in spite of including 
SN feedback. The corresponding decrease in mass in warm-hot 
gas at high densities (which can be considered as the intra-group 
and intra-cluster medium; ICM), suggests that some of the gas in 



this phase is accreted and fuels star formation. However, the exact 
evolutionary path of the gas in T - A phase space might be more 
complex than this. 

Perhaps the most remarkable result is the fact that feedback 
from AGN has a very strong impact on the thermal state of the 
diffuse gas. Indeed, comparison of the bottom panels shows that 
an additional ~ 15 per cent of the total gas mass in the simulation 
is shock-heated to temperatures above log (T / K) = 4.7 and pushed 
into regions of low density (A < 10 2 ), such that by z = 0.25 nearly 
60 per cent of the gas mass ends up in the WHIM. A comparison of 
the mass distributions among the different phases predicted by REF 
and AGN suggests that half of the additional WHIM mass, i.e. ~ 7 
per cent, at z = 0.25 is removed at higher redshifts mainly from the 
ISM (thus reducing the mass in stars at z = 0.25 by a factor ~ 3), and 
from the ICM. The remaining ~ 7 per cent of the WHIM mass ap- 
parently comes from the IGM. Comparison of the gas mass fractions 
in the warm-hot diffuse phase between the models NOZCOOL and 
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AGN suggests that SN and AGN contribute roughly a similar amount 
to the baryon content of the WHIM. Equally important, the gas mass 
in this gas phase predicted by the models NOSNJJOZCOOL (~ 30 
per cent) and AGN (~ 60 per cent) indicates that (strong) feed- 
back (both by SN and AGN) may be able to shock-heat an amount 
of gas comparable to the gas shock-heated via gravitational infall. 
This results thus indicate that it is crucial to understand feedback 
processes on super-galactic scales before being able to make any 
reliable predictions about the baryon content of warm-hot gas in the 
Universe. 

Consider finally the hydrogen neutral fraction (nni/nn) indic- 
ated by the dotted contours. The logarithmic value of (n H i/n H ) is 
indicated next to the corresponding contour in the top-left panel, and 
they are identical in all the other panels. At a fixed temperature, the 
neutral hydrogen fraction increases with density, since the ionisation 
state of the gas is dominated by photo-ionisation. However, at suffi- 
ciently high densities, i.e. A > 10 2 at z = 0.25, collisional ionisation 
dominates and the neutral hydrogen fraction depends only on the 
gas temperature, resulting in contours running parallel to the A-axis. 
In either case, the neutral hydrogen fraction steeply decreases with 
temperature at all densities (see also Richter et al. 2008 ; Danforth 
|et al.|2010") . As a consequence, the gas at densities and temperatures 
characteristic of the warm-hot diffuse gas is expected to be highly 
ionised. In particular, the model AGN predicts that the vast major- 
ity of the gas in the WHIM contains a neutral hydrogen fraction 
(nHi/nH) ;£ 10" 6 . This has important implications for the detectabil- 
ity of this gas phase via H i absorption which will be discussed in 
detail in Sec. 14*2.11 
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4.2 Observability of (warm-hot) gas using (broad) H i 
absorption 

In this section we explore to what extent the actual gas mass (distri- 
bution) in the cool and warm-hot diffuse phases are traced by the 
H i detected in absorption. In particular, we investigate the thermal 
state of the gas traced by absorbers selected in terms of their line 
width, albeit only on a statistical basis. At the same time, we invert 
the approach and inquire about the spectral signatures (and physical 
properties) of H i absorbers arising in warm-hot gas. Even though 
we are interested primarily in broad absorbers, we include narrow 
absorbers in our analysis as well. This allows for a more robust 
interpretation of our results. We define the following classes of H i 
absorbers in terms of their spectral and/or physical properties: 

(i) NLA: Hi absorber components with Doppler parameters 
i>H i < 40 km s" 1 . We adopt the notation introduced by Lehner et al. 

prjoT) . 

(ii) BLA: Hi absorber components with Doppler parameters 
£>hi > 40 kms"' that satisfy the sensitivity limit introduced by 
Richter et al. 1 2006a[l and adopted in other studies ( [Danforth et al.| 
2010[|Wiiriger et al.|2010| > 

1 >3x 10 12 (S/Nr' 



\£> H i/kms 

This limit is equivalent to a Hi Lya central optical depth 
T > 2.27 (S/NT 1 . 



(2) 



We feel that using a detection limit in terms of To is more intuitive 
than the limit in terms of (Nhi/^hi), in particular for small values 
of To, since in this limit t{v) ss 1 - F(v). Henceforth, we will use 
eq. {2} instead of the limit in terms of (N H i/£>hi) as our second BLA 
selection criterion; also, all corresponding results will be expressed 
in terms of To rather than (Nhi/^hi)- 



Figure 8. Example of a sightline through a simulation box at z = 0.25 from 
a run adopting our fiducial model (AGN). The top panel shows a synthetic 
spectrum generated by specwizard using S/N=50 (black) and the correspond- 
ing fit (red) centred on a simple BLA with £>hi ~ 63 kms" 1 , and showing 
an additional narrow H i Lya absorption line (vertical marks flag the line 
centres). The subsequent panels show, from top to bottom, the gas overdens- 
ity, the gas temperature, and the peculiar velocity along the sightline. Black, 
solid lines indicate smoothed quantities; red, dashed lines Hi optical-depth 
weighted quantities. The dotted horizontal line in the second panel (from 
the top) indicates the detection limit at S/N=50 (see eq. |A3) expressed in 
terms of A using the Nhi(A) relation discussed in Sec. |3.5.T1 The values 
in the lower-left and lower-right corner of the top panel give, respectively, 
the H i column density and the H i optical depth-weighted total hydrogen 
column density, and the overdensity and the temperature of the gas traced 
by the BLA. Note that log(7 - / K) = 5.32 corresponds to a thermal width 
bj = 59.0 kms" 1 , implying that the width of this BLA is dominated by 
thermal broadening. 



(iii) hot-BLA: BLAs with optical-depth weighted temperatures 
r H[ > 5 x 10 4 K. This class is defined in order to isolate BLAs 
genuinely tracing warm-hot gas. Note that our adopted temperat- 
ure threshold is lower than the actual temperature implied by Dop 
pier parameters b Hl > 40 km s" 1 assuming pure thermal broadening 



This is, however, not an issue since as we have shown ins Sec. |3.5 
the line width of BLAs is never entirely set by thermal broadening. 

As in previous sections, all the results presented here refer 
to simple absorbers as defined in Sec.|3.5| unless stated otherwise. 



8 The thermal width of an Hi line arising in gas at T > 5 X 10 4 K is 
b m > 30 kms" 1 . 
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Table 4. Number fraction (in percent) of simple absorbers relative to the 
total number of H i components in each class identified in synthetic spectra 
at z = 0.25 with S/N=50 for different models. 





NOSNJNOZCOOL 


NOZCOOL 


REF 


AGN 


all" 


43 


44 


44 


47 


NLA b 


37 


37 


37 


38 


BLA 


4.5 


5.5 


5.6 


7.1 


hot-BLA 


1.0 


1.7 


1.7 


3.1 



a Note that about half of the identified absorbers are single-component, 
irrespective of the model. 

b Note that the fraction of NLA is similar for all models, suggesting that 
the gas traced by these absorbers is not significantly affected by feedback. 

Tab. [4] gives an overview of the statistics of simple absorbers in 
each of the classes defined above, for all four models considered 
here. An example of such a simple (broad) Hi absorber is shown 
in Fig. [8] The top panel shows the spectrum (black) and corres- 
ponding fit (red) centred at a H i Lya- line with b Hi = 63 km s _1 and 
Nhi ~ 10' 3 cirT 2 . The next panel shows the smoothed overdensity 
(black) and the H i optical depth-weighted overdensity (red dashed) 
along the sightline. Note the relative simple density structure of the 
absorbing gas. The smoothed (black) and Hi optical depth- weighted 
(red dashed) temperatures are shown in the third panel (from the top). 
Clearly, the gas giving rise to the BLA shown in the top panel has a 
temperature T > 10 5 K. As a consequence of the high temperature, 
the absorbing gas is highly ionised, with a total H i column dens- 
ity N H ~ 4 x 10 19 cirT 2 , and it represents thus a significant baryon 
reservoir. The slight off-set between the BLA and the density peak 
of the absorbing gas is due to its (small) peculiar velocity along the 
sightline (bottom panel). 



4.2.1 Spectral sensitivity 

Before we investigate the physical and statistical properties of BLAs 
in our simulations, we need to assess how setting a fixed sensitivity 
limit as given by eq. |2]l may bias the detection of warm-hot gas. 
Under rather simple assumptions, it is possible to model the H i Lya 
central optical depth, r , of Hi absorbing gas as a function of its 
temperature and density. This allows one to put constraints on the 
physical state of the gas phase traced using H i absorption, given a 
set of instrumental limitations that lead to a minimum detection (or 
sensitivity) limit. Conversely, following this approach it is possible 
to estimate the sensitivity needed to detect gas at a given temper- 
ature and density. We describe the basic assumptions and give a 
detailed calculation of our model in Appendix |B| In particular, we 
show that it depends on the assumed size of the absorbing structure. 
With no better estimate at hand, we assume the absorbers to have a 
linear siz^Jof order the local Jeans length (Schaye 2001 , see also 
eq. |B5[ l. Note that we have already showed in Section 3.5.1| that 



this assumption can account for the Nhi(A) relation predicted by 
the simulations. Also, our model neglects non-thermal broadening, 
implying that all sensitivities in terms of t given henceforth are 
strict lower limits. 

We now investigate the relation between the gas mass distribu- 
tion in our simulations and the gas mass detected in Hi absorption. 
For each of the BLAs in the line sample obtained for each model 



Note that the simulation runs we use do resolve the Jeans length, in 
particular at the relatively low densities consider here. 



considered here, we estimate Am and T Hl , and plot the resulting 
distribution on the 7m - Am plane. The result is shown in Fig. [9] 
The coloured areas indicate the cumulative number fraction (in per 
cent) of BLAs at a given density and temperature. The gray solid 
contours correspond in each case to the gas mass distribution shown 
in Fig. [7] Note that for the contours the axes indicate the actual gas 
overdensity and gas temperature. We plot in each panel a series of 
black dashed contours which indicate the central optical depth, To, 
as a function of A H i and T H i as given by eq. pS6[ >. The corresponding 
contour values are indicated next to each curve only in the top-left 
panel, but they are identical for all the other panels. Note in particu- 
lar the thick dashed contour (magenta) which indicates our adopted 
sensitivity limit as given by eq. {2) for S/N=50, i.e. logr = -1.34 
(orlog(N H ,/Z>Hi) = 10.8). 

One notable feature in this figure is the bi-modality of the dis- 
tribution of gas traced by broad H i absorbers, irrespective of the 
model. We see in each case a population of BLAs tracing gas at 
low temperatures (7m < 5 x 10 4 K) and overdensities log A H i < 0.5, 
and a second population tracing warm-hot gas at T H i > 5 x 10 4 K 
and overdensities log Am > 0.5. Note, however, that the amplitude 
of the distribution varies from model to model. Comparing the gray 
contours to the coloured distribution we see clearly that the Hi 
detected in absorption traces only a fraction of the gas mass in 
the simulations. In particular, the model AGN (bottom-left panel) 
shows a large fraction of gas mass at 10 s K < T Hl < 3 x 10 5 K and 
log Ah i ~ 0.5 which is not detected in Hi absorption. The same is 
true for the models NOZCOOL and REF, although at slightly dif- 
ferent temperatures and overdensities. Consideration of the thick 
dashed contour reveals that this is a selection effect, i.e. the gas is 
simply not detectable at our adopted sensitivity limit. As discussed 
above, this comes about because the gas at such high temperatures 
and relative low densities is extremely ionised, with neutral hy- 
drogen fractions log (n H i/n H ) ^ -6, and its absorption simply falls 
below our adopted detection threshold (cf. Fig.jBTJ. 

Thus, in our fiducial model BLAs detected in spectra with 
S/N=50 typically have -2 < log t < 0, while the bulk of the gas 
mass at 10 5 K < T m < 3 x 10 5 K and log A H i ~ 0.5 is predicted to 
give rise to absorption with log To < -2. The fact that we do detect 
in absorption some of the gas at temperatures and densities which 
correspond to sensitivities slightly smaller than our adopted value 
(i.e. to the left of the thick dashed contour) reflects the simplicity of 
the assumptions that go into modelling the absorption strength in 
terms of A and T. Nevertheless, the expected and actual detections 
are fairly consistent with each other. Based on this, we estimate 
that in order to detect most of the baryonic mass in the WHIM 
using thermally broadened H i absorption, spectra with very high 

S/N are required that allow detection at the log To 2 level, which 

is roughly an order of magnitude lower than the typical sensitivities 
adopted in BLA studies ( |Richter et al.|2006a||Danforth et al.|2010| 
|Willigeretal.|2010) . 

4.3 BLA number density 

Since our adopted sensitivity limit matches the value used to identify 
BLA candidates in real QSO spectra, we may directly compare the 
predicted and observed line frequencies. Applying the selection cri- 
teria described above (i.e. &m > 40kms~' and To > 2.27 (S/N) -1 ) 
to our AGN H i line sample obtained from spectra with S/N=50 at 
z = 0.25 results in 6120 BLA candidates, which corresponds to a 
line-number density (dA7dz) BLA = 33 ± 6, where the quoted uncer- 
tainty is pure Poissonian. The number densities of BLAs identified 
in spectra with S/N=30 and S/N=10 are given in Tab. [5] For com- 
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Figure 9. Distribution of temperature and overdensities of the gas traced by simple BLAs (coloured regions) identified in spectra with S/N=50 at z = 0.25 
obtained from various models. In each panel, the coloured regions show the cumulative number fraction (in per cent) of BLAs. For instance, the red region 
contains 25 per cent of the total number of BLAs, while the green and red regions together contain half of them. The gray contours correspond to the gas mass 
distributions shown in Fig. [7] Note the change in scale on both the x- and y-axes with respect to Fig. [7] The vertical (horizontal) solid line in each panel indicates 
the density (temperature) threshold adopted to distinguish the different gas phases (see Tab.|3l. The dashed contours, which are identical in all panels, indicate 
the central optical depth, to, as a function of (over-)density and temperature (see Appendix[B] eq. |B6( . For clarity, the corresponding logarithmic value of To is 
included next the each curve only in the top-left panel. The magenta dashed contour indicates our adopted sensitivity limit for S/N=50, log To = - 1 .34 (eq.[2j. 
The percentage (orange) in each panel indicates the fraction of baryonic mass of the gas traced by simple BLAs with 7"h i > 5 x 10 4 K, i.e. simple hot-BLAs. 
Reassuringly, the overwhelmingly majority of the identified BLAs lie below the thick magenta contour, irrespective of the model. However, the large fraction of 

mass in gas at T > 10 5 K and A < 10 (in particular for AGN) is not observable at our adopted sensitivity. The dashed contours indicate that at least log To 2 

is required. BLAs selected in terms of their width and eq. j2j thus only trace the low-temperature regime of the warm-hot diffuse gas, since the gas at higher 
temperatures is not observable, probably due it its high ionisation degree (see Fig. |Bl[ . Note that the estimated baryon content of the absorbing gas in each case 
corresponds to ~ 10 per cent of the actual baryonic mass in the this phase for the corresponding model (cf. Fig.|JJ. 



pleteness, we also include in this table the corresponding numbers 
for NLAs. Note that the absolute number of BLAs decreases while 
the relative number of simple BLAs increases with decreasing S/N. 
This is due to a low S/N causing line blending and increasing the 
line widths, as discussed in Sec. |3.3| (see also R ichter et al.| 2006a 
their figure 2). 

In Fig. [To] we compare the cumulative BLA number density 
as a function of r predicted by our fiducial run obtained from 
spectra with various S/N values to the available observational results 
obtained from QSO spectra at comparable redshifts and with similar 
(average) S/N. The arrows indicate the BLA sensitivity limit in 
terms of To for the corresponding S/N value as given by eq. {2). The 
predictions for simple BLAs are indicated by the solid lines. Since 
our definition of simple absorbers is somewhat arbitrary, we include 
the corresponding predictions for all, i.e. simple and complex, BLA 



candidates as well (dashed lines). These two sets of lines thus span 
our predicted line-frequency range for S/N € [10, 50]. 

It is noteworthy that our predictions are broadly consist- 
ent with the observed range of BLA number densities. For 
example, our prediction for S/N=30 (blue solid) agrees with 
the results by |Lehner et al.| j2007| green data points). These 
authors find in their data with an average S/N as 15, discard- 
ing lines with associated line-parameter errors larger than 40 
per cent, a fraction of single-component BLAs close to 30 
per cent and a mean (averaged over 7 sightlines) line-number 
density (dA'/dz) BLA = 30 ± 7 at log t > -1.10 (equivalent to 
logOW&m) S 1L02), and (cW/dz) BLA = 9 ± 3 at logT > -0.66 
(or log (N m /b m ) > 11.46 ). Both our S/N=50 (black) and 
S/N=30 (blue) predictions match the result by Danfort h et al.| 
l |2010| orange data point), who find (d/V/dz) BLA = 18 ± 11 at 
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Table 5. Line-number density, (dN/dz), and corresponding Poisson uncer- 
tainties for BLAs and NLAs identified in spectra with various S/N values 
from our fiducial run at z = 0.25. 





S/N=50 


S/N=30 


S/N=10 


all" 


460±21 


332±18 


144±12 


simple 


214±15 


175±13 


94±10 


NLA (simple) 


173+13 


140±12 


73±9 b 


NLA (all) 


343±19 


248±16 


110+10 b 


BLA (simple) 


33±6 


28±5 


15±4 


BLA (all) 


95±10 


66±8 


25 ±5 



a This corresponds simple and complex H i absorbers taken together. 
b For reference, |Lehner et al.|{2007| find a mean (dAVdz)NLA = 66 ± 17 
over 7 sightlines for all NLAs in their data with an average S/N » 10, 
discarding lines with associated Voigt-profile parameter errors larger than 40 
per cent. 



(Nhi, £>hi) using our fiducial model. We now focus on the physical 
conditions of the gas giving rise to H i absorbers identified as simple 
BLAs in our synthetic spectra at z = 0.25; these correspond to single- 
component H i absorbers falling within the region defined by the 
polygon in Fig. [6] which is defined through two criteria: a) the line 
width satisfies £> H i > 40kms _I ; b) the central optical depth obeys 
eq. l(2j, i.e. T > 2.27 (S/N) -1 . 

Given the importance of broad Hi absorbers as potential 
WHIM tracers, and the dependence of the predicted WHIM mass 
fraction on the adopted physical model, we next explore the rela- 
tion between the measured line width and the temperature of the 
absorbing gas using the models introduced in Sec. [2] We further use 
our fiducial model (AGN) to investigate the ionisation state, neutral 
hydrogen fraction, and metallicity of the gas traced by different 
classes of H i absorbers. Finally, we estimate the baryon content of 
the gas traced by BLAs using all models, but compare our results to 
observations using only the model AGN. 



AGN;S/N=50 

^ l-» H> S/N=30 

S/N=10 

'*\ simple 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

logx 

Figure 10. Cumulative number density of BLAs as a function of central 
optical depth, To. Measurements available in the literature are plotted along 
with the predictions of our fiducial run AGN at z, = 0.25 (see text for details). 
The solid lines show the result for simple BLAs identified in 5000 spectra 
with S/N=50 (black), S/N=30 (blue), and S/N=10 (red). We include our 
results for all (i.e. simple and complex) BLA candidates as well (dashed 
lines). The arrows indicate the BLA sensitivity limit in terms of To for the 
corresponding S/N value as given by eq. j2j. 

log i"o ^ -1.12 (log (Nhi/^hi) > 11.0) along seven sightlines at 
z < 0.5, spanning a total redshift path Az = 2.193, with S/N > 5. 
Their sample consists of 15 single-component BLAs and 48 
BLAs with uncertain velocity structure. Wil liger et al.| ( |2010| 
magenta data point) estimate their detection limit at log t > - 1 .22 
(or log(N H ,//j H i) £ 10.9) and obtain (d/V/dz) BLA = 25+15 and 
(d/V/dz) BLA = 8.5 ± 8.5 for their full and single-component BLA 
samples, respectively, along a single sightline (Az = 0.117) and 
corresponding spectrum with S/N=20 - 30. Finally, Richt eFet al.| 
(2006a, black data point) measure (dA'/dz) BLA = 22 using their 
reliable sample of single-component BLAs, detected along four 
sightlines at z < 0.4 (Az = 0.928), at a sensitivity log t > -0.82 (or 
log (Nhi//?hi) ^ H-3), in spectra with an average S/N « 15. 

4.4 Physical properties of broad H i absorbers 

In Sec.[3]we discussed the relation between the physical conditions 
of the gas traced by typical H i absorbers and their line observables 



4.4.1 Temperature distribution of (broad) lines 

In Fig.[TT]we show the distribution of temperatures of the gas traced 
by both BLAs (red histograms) and NLAs (blue histograms) identi- 
fied our spectra with S/N = 50 at z = 0.25 obtained from different 
models. The dashed, vertical line in each panel indicates the tem- 
perature threshold adopted to separate cool from warm-hot gas (see 
Tab. [3). 

In all models, the temperature distribution of the gas traced 
by NLAs shows that lines with Doppler parameters ini < 40kms~' 
preferentially arise in gas at temperatures T H i S5x 10 4 K with a 
peak at T H i ~ 10 4 K, as expected from their width. The few lines 
which are narrower than allowed by the temperature of the absorbing 
gas (i.e. the section of the blue histogram to the right of the ver- 
tical, dashed line) are mostly weak lines that fall below the formal 
completeness limit, as discussed in Sec. |3.5,2| But some of these 
lines are real detections, which suggest that gas at different temper- 
atures overlaps in velocity space (due to redshift-space distortions), 
and some of it may even overlap in position space, indicating the 
existence of multi-phase absorbing structures. In models with feed- 
back (NOZCOOL, REF, AGN), the temperature distribution of the 
gas traced by NLAs is broader than in the model without feedback 
(NOSN_NOZCOOL), and the fraction of lines with 'un-physical' 
widths is higher. This can be explained as follows. Outflows driven 
by SNe and AGN follow the path of least resistance in space, thus 
escaping into the voids while leaving the cooler, denser filaments 
intact < |Theuns et al.|2002| >. With increasing feedback strength, the 
cross-section of such high-temperature outflows increases as well, 
and so does the chance for a random sightline to intersect both cool, 
dense filaments and shock-heated material, with their corresponding 
absorption overlapping from time to time along the spectrum. 

Note that the NLA temperature distributions in the models 
NOZCOOL and REF are very similar to each other, and the same is 
true for the corresponding BLA temperature distributions. Moreover, 
the fraction of hot-BLAs is only slightly lower in REF than in 
NOZCOOL. This indicates that metal-line cooling is of secondary 
importance in setting the thermal state of the gas phase traced by 
(broad) Hi absorbers. 

In contrast to narrow Hi absorbers, broad Hi lines trace gas 
around two different temperatures, r H[ ~ 10 4 K and T Hl ~ 10 5 K, 
irrespective of the model (red histograms). Clearly, BLAs arising in 
gas at low temperatures must be subject to substantial non-thermal 
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Figure 11. Temperature distribution of the gas traced by simple Hi absorbers identified in spectra at z = 0.25 with S/N=50 obtained from different models. The 
blue histograms show in each case the distribution for gas traced by narrow (i.e. &hi < 40kms~') absorbers, while the red histograms indicate the temperature 
distribution of gas traced by BLAs. The percentage shown in each panel indicates the number fraction of BLAs arising in gas at 7hi ^ 5 X 10 4 K, i.e. hot-BLAs. 
Note the correspondence between increasing mass fraction of warm-hot diffuse gas in Fig.|7]and the increasing fraction of BLAs tracing warm-hot gas shown 
here. The distribution predicted by our fiducial run (AGN) suggests that, statistically, two out of three BLA candidates trace gas at T > 5 X 10 4 K. 



broadening, such as bulk flows and/or Hubble broadening] Since 
the line width of an Hi absorber with £> H i = 40km s arising in 
gas at T ~ 10 4 K is completely dominated by non-thermal broad- 
ening, its linear size (assuming the line width is entirely due to 
Hubble broadening) must be ~ 500 kpc, which is consistent with 
the Jeans length of a filament with a mean density nn 
and temperature T ~ 10 4 K (eq. 



B5i 



Interestingly, the fraction of BLAs tracing gas at low 
(r H i ~ 10 4 K) and high (r H i ~ 10 5 K) temperatures is very different 
in each model, as can be judged qualitatively by the shape of the 
corresponding histograms and quantitatively by the percentage in- 
cluded in each panel which gives the fraction of hot-BLAs, i.e. BLAs 
tracing gas with T > 5 X 10 4 K. In fact, the ratio of BLAs tracing 
cool gas to those tracing warm-hot gas appears to be very sensitive 
to feedback strength. In principle, this could be used to constrain 



Note that our simulations lack the resolution to capture small-scale turbu- 



lence within the gas. 



feedback models observationally. The caveat is that a statistically 
significant sample of confirmed BLAs would be required for which 
the gas temperature can be measured reliably. In our fiducial model 
(AGN), which includes the strongest feedback, the majority of the 
BLAs trace gas at T Hi ~ 10 5 K, with little contamination by non- 
thermally broadened lines. In fact, two out of three BLA candidates 
arise in gas at temperatures T > 5 x 10 4 K. 

The bi-modal character of the gas temperature distributions for 
BLAs predicted by the model NOSNJJOZCOOL is also consistent 
with previous results. Ric hter et al.| ( [2006b[ l find in a simulation 
that included a model similar to our NOSNJJOZCOOL that ~ 30 
per cent of the BLAs trace gas at T < 1 x 10 4 K, and a significant 
fraction gas at 7hi > 5 x 10 4 K. The quantitative difference between 
theirs and our result for NOSN _NOZCOOL is probably caused by 
the difference in the simulation methods used. Using a linear model 
to fit the thermal to total line width, |Richter et al.|p006b| > find that 
br/byn = 0.91 (no error quoted). We note that we do not find such 
a tight correlation between b T and bm, but if we perform a linear 
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Figure 12. Distribution of various physical quantities characterising the gas giving rise to single-component absorption features identified as NLA (blue), BLA 
(red), and hot-BLA (orange) in spectra at z = 0.25 with S/N=50 obtained from our AGN model. Left: Neutral hydrogen fraction Middle: Total hydrogen column 
density. Right: Metallicity relative to the solar value (Z Q = 0.0127). The vertical dashed lines in each panel indicate the median values of the corresponding 
distributions. On average, NLAs and (hot-)BLAs appear to trace different gas phases. 



fit between these two quantities, we fincplfer/fem = 0.641 ± 0.004. 
This is consistent with the result presented in Sec. |3.5.2| that thermal 
broadening on average contributes with (at least) 60 per cent to the 
total line width of BLAs. 

|Lehner eTaf] ( [2007] > argue that broad Hi lines may trace 
both cool and warm-hot gas, but that the majority of BLAs trace 
gas at T ~ 10 5 - 10 6 K if their width is dominated by thermal 
broadening. As mentioned in the previous paragraph, thermal 
broadening accounts for a significant fraction to the total line width 
of single-component, broad Hi absorbers. Thus, our simulations 
are consistent with the result inferred from observations that these 
absorbers do preferentially trace gas at high temperatures, at least in 
models with (some type of) feedback. 

In summary, we find that in the absence of feedback BLA 
samples are contaminated by a large fraction of non-thermally 
broadened lines. Conversely, the fraction of broad H i absorption 
lines tracing gas at temperatures 7hi ~ 10 5 K increases when feed- 
back is included. For instance, our fiducial model predicts that, in 
a statistical significant sample, 67 per cent of the BLAs trace gas 
at 7hi > 5 x 10 4 K. Our results thus strongly support the idea that 
reliable BLAs detected in real absorption spectra are genuine tracers 
of gas at such high temperatures. 

4.4.2 Neutral hydrogen fraction, total hydrogen column density, 
and metallicity 

In Fig. [12] we show the distribution of neutral hydrogen fraction 
(left panel), total hydrogen column density (middle panel), and 
(local) metallicity (right panel) of the gas traced by narrow absorbers 
(blue), BLAs (red), and hot-BLAs (orange) identified in spectra with 
S/N=50 obtained from our fiducial model. The vertical dashed line 
indicates in each case the corresponding median value. 

In general terms, the physical properties of the gas traced by 
BLAs and hot-BLAs show similar distributions and comparable me- 
dian values, but they are somewhat different from the corresponding 
properties of the gas traced by NLAs. For example, the median neut- 
ral hydrogen fraction of gas traced by NLAs is log(nm/nii) 5.5, 

which is slightly higher than the neutral hydrogen fraction of the 

gas traced by (hot-)BLAs, log(n H i/n H ) 6. This is expected since, 

as we have shown previously, the temperature of gas giving rise to 

For reference, the corresponding results for S/N=30 and S/N=10 are 
h T /hn, = 0.622 ± 0.04 and b T /b Hl = 0.556 ± 0.006, respectively. 



broad H i absorption is, on average, higher than the temperature of 
gas giving rise to narrow H i absorbers. Furthermore, the median 
total hydrogen column density of the gas detected via (hot-)BLAs 
is N H ~ 6 x 10 18 cm~ 2 , which is several times larger than the me- 
dian total hydrogen column density of the gas traced by NLAs, 
Nh ~ 2 x 10 18 citT 2 , and its distribution extends out to significantly 
larger values, N H ~ 10 20 cm~ 2 , as compared to N H ~ 10 19 cnr 2 . This 
is due to several factors. First, as shown in Figs. |Bl| and |9[ high- 
temperature gas detected via H i absorption at a fixed sensitivity 
necessarily has a higher density with respect to gas at lower tem- 
peratures. Also, a higher density implies an average higher Nm as 
a consequence of the Am - Nm correlation. Finally, gas at high 
temperature has a lower neutral hydrogen fraction, which in turn 
yields higher total hydrogen column densities for a given Nm- The 
high(-er) total hydrogen density of the gas traced by broad H i ab- 
sorbers implies that its baryon content is considerable. We will come 
back to this point in more detail in Sec. |4.5| 

Quite interesting is the difference between the gas metallicity 
distributions. While the metallicity of the gas traced by NLAs shows 
a broad distribution with a tail extending to very low values and a 
median Z ~ 10~ 3 Z o , the metallicity distribution of the gas traced 
by broad absorbers is narrow, with most values falling in the range 
(0.001, 1) Z , centred around Z ~ 10~ 2 Z o . On average, the metalli- 
city of the gas traced by (hot-)BLAs exceeds the metallicity of the 
gas traced by NLAs by an order of magnitude. 

These results together indicate that broad H i absorbers trace 
gas that is physical distinct from the gas traced by narrow H i ab- 
sorbers, as already mentioned in Sec. [333] In particular, the rel- 
atively high level of enrichment is inconsistent with the idea that 
BLAs trace primordial gas that is sinking along filaments towards 
the centre of high-density regions, as commonly assumed. Rather, 
our results suggest that broad H i absorbers may be tracing recent 
(or on-going) galactic outflows, and/or gravitationally shock-heated 
gas that has been enriched by galactic ejecta at early epochs. 

4.5 Baryon content of H i absorbing gas 

In this section we briefly investigate the dependency of the predicted 
baryon fraction of the gas traced by H i absorbers on the adopted 
physical model. As we have done in Paper I for O vi absorbers, we 
estimate the baryon fraction, i.e. the baryon density relative to the 
critical density p c , in Hi absorbers using 

( NLOS \~' "LOS Wabs N 

n b (Hi)=— ttZ a *; EE t^t^^V^' (3) 

Pc { H o ) jrf jbi (X »)hi (n H i/n H )Hi 
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Figure 13. Top: Baryon content (in per cent) relative to the cosmic value 
fib = 0.0418 of the gas traced by H i absorbers identified in synthetic spectra 
with S/N=50 for different model runs. The H i absorber sample has been 
dissected into NLA (blue squares), BLAs (red filled circles), and hot-BLAs 
(orange open circles); these classes have in turn been divided into single- 
component absorbers (solid lines) and complex systems (dotted lines). The 
black, dot-dashed lines indicate in each case the result for all (i.e. simple and 
complex) absorbers of a given class. Bottom: Gas mass traced by a given 
Hi absorber type relative to the actual gas mass in the phase expected to 
be traced by that particular type (see text for details). Symbols, lines, and 
colours as in the top panel. 



where m H is the hydrogen mass, and X H and (n H i/n H )Hi are the 
optical-depth weighted hydrogen mass fraction and neutral hydrogen 
fraction, respectively. Note that Nm, (Xh)hh an d (nm/nH)Hi are 
computed for individual absorbing components along each sightline, 
but we have omitted the running indices for simplicity. 

The top panel of Fig. [TJJ shows the baryonic mass fractions, 
£lb(Hi)/£l b , predicted by various models in different types of ab- 
sorbers: NLAs (blue squares), BLAs (red filled circles), and hot- 
BLAs (orange open circles), where each of these classes have been 
sub-divided into simple (solid lines) and complex (dotted lines) 
absorbers (see Sees. |3.5] and |4.2| . Note that we consider the ba- 
ryonic mass fraction of the gas traced by both simple and complex 
absorbers, since our adopted criterium to define simple absorbers 
is somehow arbitrary. The net baryon fractions of simple and com- 
plex absorbers (of a given class) taken together are indicated by the 
dot-dashed lines. 

The baryonic mass fraction in NLAs (simple or complex) is 
very similar in all models, being only slightly lower in our fiducial 



run AGN. This is consistent with the lower gas mass fraction in 
the cool diffuse gas (which is expected to be traced by NLAs) in 
this model compared to all other models (see Fig. [7jl. The small 
difference in £l b (NLA)/n b between all the models indicates that 
feedback has a negligible impact on the gas phase typically traced 
by NLAs. This in turn is consistent with the fact that the predicted 
H i statistics (which are dominated by these absorbers) are rather 
insensitive to variations in the feedback model (see Appendix O, 

The baryonic mass fraction of gas traced by simple BLAs is 
relatively low, and varies significantly between the models, from ~ 3 
per cent {NOSN JNOZCOOL) to ~ 7 per cent (AGN). The baryonic 
mass fractions in complex BLAs are much higher than the baryonic 
mass fractions in simple BLAs, and they are also very different in 
each model. For instance, the baryon fraction of complex BLAs in 
the model NOZCOOL is higher (~ 5 per cent) than in the model 
NOSNJJOZCOOL; this is consistent with the fact that SN feedback 
significantly increases the mass in warm-hot diffuse gas, as has been 
shown previously (see Fig.|7]and corresponding text), together with 
the idea that BLAs preferentially trace this gas phase. In contrast, the 
baryon fraction in complex BLAs predicted by the models REF and 
AGN, is lower (by ~ 5 and ~ 10 per cent, respectively) compared to 
predictions of the model NOZCOOL. 

The lower baryon fraction in complex BLAs predicted by the 
model REF with respect to NOZCOOL can be explained as follows. 
Complex absorbers trace kinematically disturbed gas, most probably 
SN-driven outflows. These ejecta carry heavy elements with them, 
which allow a significant fraction of the gas to cool down radiatively, 
thus reducing the number of thermally broadened lines and their 
net baryonic mass. However, the lower baryonic mass content of 
complex BLAs in the model AGN with respect to all other models 
is in contrast with the actual total mass fraction in the warm-hot 
phase predicted by this model, which is higher compared to all other 
models (see bottom-left panel of Fig.[7]l. The discrepancy between 
the predicted mass fraction in the warm-hot phase and the baryon 
content of the gas traced by BLAs in the AGN model can be under- 
stood as consequence of the limited sensitivity. As already discussed, 
AGN feedback shifts a significant fraction of gas into the warm- 
hot phase; however, most of this mass ends up at temperatures and 
densities which lead to a H i fraction and corresponding absorption 
signal that is beyond our adopted detection limit (see Fig.[9]l. 

Note that the baryon fractions traced by hot-BLAs are only 
sightly lower than for BLAs, irrespective of the model. This is im- 
portant because it implies that the contamination of the BLA sample 
by non-thermally broadened lines does not significantly affect the 
inferred baryon fraction of the WHIM. In other words, the baryon 
fraction in warm-hot gas is dominated by the absorbers arising in 
gas at the highest temperatures. This is a direct consequence of the 
steep decline of the hydrogen neutral fraction with temperature. 

The bottom panel of Fig.fTJJshows the recovered mass fraction 
(in per cent) of a given gas phase. This quantity is defined as the 
total baryonic mass in a given absorber class relative to the actual 
gas mass in the phase (Fig. [7) expected to be traced by that particular 
absorber class. So, for example, the gas mass recovered from hot- 
BLAs is given in each model as the total baryonic mass in these 
absorbers (orange, filled squares in the top panel of Fig. 1 1 3|> divided 
by the gas mass in the warm-hot diffuse phase (orange percentages 
in the top-left sections of Fig. |7J. The gas mass recovered from 
NLAs is correspondingly given as the total baryonic mass in these 
absorbers (blue squares in the top panel of Fig. 1 1 3| > divided by the 
gas mass in the cool diffuse phase (percentages indicated in the 
bottom-left sections of Fig.|7j(. 

Apparently, none of the absorber classes traces the total gas 
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mass in its corresponding (expected) phase, with exception per- 
haps of the full (i.e. simple and complex) BLA sample (red filled 
circles) in the model NOSN_NOZCOOL. The gas mass fraction 
traced by both simple and complex NLAs is very similar in all mod- 
els, and taken together these absorbers trace between ~ 70 per cent 
(NOSN_NOZCOOL) and ~ 90 per cent (AGN) of the gas mass in 
cool diffuse gas. This again is consistent with our previous state- 
ment that the gas in this phase is left almost intact by feedback 
mechanisms such as SN-driven winds and AGN outflows, 

Note that simple (hot-)BLAs trace roughly 10 to 15 per cent of 
the true baryonic mass in warm-hot gas, irrespective of the model. 
This suggests that baryonic mass estimates based on this type of 
absorbers are robust. In contrast, the gas mass traced by complex 
(hot-)BLAs is very different in each model. As already mentioned 
above, in the model NOSN-NOZCOOL the full BLA sample traces 
practically all the mass contained in warm-hot gas, with complex 
BLAs contributing with more than 80 per cent to the recovered 
gas mass. This suggests that the bulk of gas shock-heated by 
gravity (which is the only possible heating mechanism in the model 
NOSNJJOZCOOL) can be fully accounted for using (simple and 
complex) BLAs, at our adopted sensitivity. The recovered WHIM 
mass is, however, systematically lower in models that include SN 
and AGN feedback, and metal-line cooling. 

Taking the result from our fiducial run AGN at face value, we 
estimate the total baryon content in gas traced by H i in our simula- 
tion at z = 0.25 to be n b (Hi)/fi b = 0.57 (S/N=50), 0.48 (S/N=30), 
and 0.29 (S/N=10). The last two values are in remarkable agreement 
with the results from observations at comparable sensitivity. Assum- 
ing a simple ionisation model, Penton, Stocke & Shull (2004) meas- 
ur^]n b (Hi)/n b = 0.31 ± 0.04 at z ~ for absorbers with column 
densities 12.5 < log (N Hl /cirT 2 ) < 17.5 andb H i < 100 kms~'. Sim- 
ilarly, a ssum ing the gas to be isothermal and photo-ionised, |Lehner| 



et al.| ( |2007| l obtain^] n(Hi)/n b =0.40 from their data with 
an average S/N « 15 and for 12.4 < log (n H i/citT 2 ) < 16.5 and 
bm < 150kms-'. 



4.5.1 Baryon content of warm-hot gas at low z 

Estimates of baryonic mass contained in the WHIM based on broad 
H i absorbers detected in real QSO spectra are very uncertain, even 
with a reliable sample of BLA candidates at hand, since they are 
highly sensitive to the ionisation state of the absorbing gas (see 
eq. [3j, which is probably dominated by collisions between ions and 
electrons in the plasma. In this case, the neutral hydrogen fraction 
is a steeply decreasing function of temperature, and an accurate 
estimate of the WHIM baryon content thus relies on a precise meas- 
urement of the temperature of the absorbing gas. As we have shown 
in Sec. |3. 5.2"] temperature estimates from the line width of broad H i 
absorbers may yield values that are uncertain by, at least, factors of 
a few. 

The first attempt to measure the baryon content of the 
WHIM using BLAs was undertaken by Ric hter et al.| (2004), 
who found f2 b (BLA) < 3.2 x 10~ 3 (A/0.73) -1 assuming CIE, 
which represents less than 8 per cent of the cosmic baryon 
budget. In a follow-up study, Richt er et al.| ( |2006a[ > measured 



£l b (BLA) > 2.6 x 10~ 3 (h/0.73)' 1 , corresponding to at least 6 per 
cent of the baryons in the Universe. These authors also assumed 
CIE, but recognised the potential importance of photo-ionisation 
(PI) in determining the ionisation state of the WHIM, and concluded 
that their baryon content measurement could be underestimated by 
15-50 per cent. 

Using a significantly larger sample than previous studies, 
|Lehneretal.| l (2007) report n b (BLA)/fi b = 0.08 assuming the gas to 
be in collisional ionisation equilibrium (CIE) for absorbers with 
13.2 < log(N H ,/cm- 2 ) < 16.5 and 40kms-' < bm < 150 kms" 1 . 
Using the same sample and assuming a hybrid model (including 
photo- and collisional ionisation) to compute the neutral hydrogen 
fraction, these authors find £l b (BLA)/n b = 0.20. Both estimates are 
based on a series of assumptions. First, in order to account for the 
possible contamination of their sample with lines broadened by unre- 
solved velocity structure or any other non-thermal mechanism, these 
authors randomly discard one third of the BLAs in their sample. 
Moreover, they assume the thermal width to be 90 per cent of the 
observed line width, based on the results from previous simulations 
by |Richter et al.| ( [2006b[ l. If, instead, the line width is dominated 
by thermal broadening, they get fl b (BLA)/f2 b = 0.13 (CIE) and 
£l b (BLA)/n b = 0.32 (PI+CIE). 

In a more recent studjFj 
port n b (BLA) = 6.0Zol " ' n ^ 
£l b (BLA)/n b 



Danforth 



et al. 



1 2010 1 



re- 
to 



14 +003 



0V0.73) , equivalent 
These authors analyse in detail the 
systematic uncertainties that afflict their (and others') baryon 
estimates, such as unresolved velocity structure, sample com- 
pleteness, ionisation corrections, and the assumed relation 
between line width and gas temperature, and find their es- 
timate to vary between fi b (BLA) = 2.3 x 10~ 3 (A/0.73)" 1 and 
£l b (BLA) = 15.2 x 10~ 3 (A/0.73)" 1 , i.e. between ~ 6 and ~ 36 per 
cent of the cosmic baryon budget. Clearly, there is still a high 
uncertainty in the estimate of the baryonic mass traced by observed 
BLAs. 

The results from our fiducial model are as follows. If we take 
simple and complex BLAs together, we find that they trace ~ 25 
(S/N=50), ~ 20 (S/N=30), and ~ 10 (S/N=10) per cent of the total 
baryon budget in our simulation. For comparison, BLAs (simple and 
complex) that arise in gas at T > 5 x 10 4 K yield ~ 24 (S/N=50), 
~ 18 (S/N=30), and ~ 9 (S/N=10) per cent, which are very close 
the values obtained from the whole BLA sample. If we restrict the 
BLAs to be single-component, the resulting baryon fractions in 
these absorbers are ~ 7 (S/N=50), ~ 6 (S/N=30), and ~ 5 (S/N=10) 
per cent. This confirms that contamination of the BLA sample by 
non-fhermally broadened lines does not significantly affect the 
inferred baryon fraction of the WHIM. Incidentally, this suggests 
that there should be little overlap between the estimates of £l b (NLA) 
and n b (BLA). 

One final important remark. We have demonstrated that the 
broad Hi absorbers trace only a fraction of the total mass in the 
WHIM phase. Thus, even in the case that one could accurately 
estimate the baryon fraction in a given sample of absorbers, there is 
still a large gap between the observed and true mass contained in 
this gas phase. Although our simulation suggests that the baryonic 
masses estimated from observations represent 1/10 to 1/3 of the true 



12 Relative to £l b = 0.0418, rather than Cl b = 0.047 assumed by |Penton etaf] 
(2004) , and re-scaled to h = 0.73. 

' 3 Relative to £l b = 0.0418, rather than £2 b = 0.044, and re-scaled to 
h = 0.73 rather than h = 0.7. 



We note that there is an error in the computation of the total absorption 
path length in Danforth et al. (2010 their Table 1), which might have 
affected their reported estimate of the baryon content in BLAs (C. Danforth, 
private communication). 
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baryonic mass in the WHIM, it is not clear how to use the measured 
baryonic mass to infer the true, total amount of baryons in this gas 
phase. 



5 SUMMARY 

In this paper, we have used a set of cosmological simulations 
from the Overwhelmingly Large Simulations (OWLS) project 
(Schaye et al. 2010) to study the physical conditions of the gas 
traced by Broad H i-Lya Absorbers (BLAs) with low and moderate 
column densities (log (Nhi/chT 2 ) < 15) observed in QSO spectra. 
We have chosen the AGN model of the OWLS suite to test the 
predictions of our simulations against a set of H i observables. We 
have investigated the impact of metal-line cooling, kinetic feedback 
by supernovae (SNe) explosions, and feedback by active galactic 
nuclei (AGN) on the distribution of the gas mass over different 
phases such as the photo-ionised intergalactic medium (IGM) and 
the shock-heated, warm-hot intergalactic medium (WHIM). Finally, 
we have explored the relation between the physical state and the 
baryon content of these gas phases and both narrow H i absorbers 
(NLAs) and BLAs. 

Our detailed results can be summarised as follows: 

• Accretion shocks due to gravitational infall into the potential 
wells of dark matter halos heat ~ 30 per cent of the total gas mass 
to temperatures r^5xl0 4 Kbyz; 



4.1 



: 0.25 (Sec. 

• Feedback by SNe and AGN each remove a similar amount of 
gas from the ISM in haloes at early epochs and displace it to the 
warm-hot diffuse phase, increasing its total mass fraction by another 
~ 30 per cent by z = 0.25 (Sec. |4,l} ; in other words, roughly half of 
of the gas mass predicted to be in the WHIM at low redshift (~ 60 
per cent) has been heated by accretion shocks while the other half is 
due to strong feedback. 

• The predictions from our simulations are in excellent agree- 
ment with standard H i observables (CDDF, line-width distribution, 
bm _ N H i correlation; Secs. |3.2| - |3,4[ >; these observables are rather 
insensitive to feedback and/or metal-line cooling (Appendix |C). 

• The line-number density of narrow (£> H i S 5 40kms _1 ; NLA) 
and broad (£>hi 3 s 40 km s ; BLA) Lya- absorbers predicted by our 
fiducial run AGN are in broad agreement with the corresponding 
observed line-frequencies (Sec. |4.3[ l. 

• The density of the H i absorbing gas shows a tight correlation 
with the H i column density, which agrees well with the analytic 
prediction of Schaye (2001 . Sec. |3. 5.1) ; this implies that our simula- 
tions are consistent with the assumption that typical H i absorbers 
are self-gravitating clouds in hydrostatic equilibrium with linear 
sizes of the order of the local Jeans length. 

• The temperature of the H i absorbing gas correlates well with 
the Hi line width for log (N H] /cm _2 j > 13, but it is a poor indicator 
of the thermal state of the gas for lower column densities; thermal 
broadening contributes, on average, with at least 60 per cent to the 
line width of BLAs (Sec. [331) . 

• The overwhelming majority of NLAs are found to trace gas at 
T ~ 10 4 K; their number, temperature distribution and baryon con- 
tent is very similar in models with/without feedback, thus strongly 
suggesting that feedback has a negligible impact on the cool, diffuse 
gas (i.e. the IGM; Sec. [44). 

• BLAs trace gas both at T ~ 10 4 K and at T ~ 10 5 K; our fidu- 
cial model, which includes feedback by both SNe and AGN, predicts 
that 2 out of 3 BLAs arise in gas at T > 5 x 10 4 K; the number ratio 



of thermally to non-thermally broadened H i absorbers is very sens- 
itive to (the adopted) feedback (model), and could in principle be 
used as an indicator of feedback strength (Secs. |4~2]|4.4} . 

• The ionisation state, the total hydrogen content, and the level of 
enrichment of the gas traced by BLAs indicates that these absorbers 
arise in gas that is physically distinct from the gas traced by NLAs; 
we argue that BLAs mostly trace gas that has been recently shock- 
heated and enriched by outflows (Secs |3.53]|4.4[ l. 

• While models including SN and AGN feedback predict a higher 
fraction of gas mass to be in the warm-hot diffuse phase, the baryon 
fraction of the gas inferred from BLAs in these models is lower 
compared to a model without feedback; the reason is that much of 
the mass is displaced to temperatures and densities for which the 
H i fraction is too low for the gas to be detectable at the adopted 
sensitivity (Sec, |4.5| >. 

• The baryon fraction of the gas traced by both NLAs and BLAs 
predicted by our fiducial model shows broad agreement with corres- 
ponding measurements from observations (Sec. |4,5[ >. 

• Baryonic mass estimates using simple BLAs are robust; in 
contrast, the gas mass traced by complex BLAs is very sensitive to 
the adopted (feedback) model (Sec. [43). 

• Our fiducial model predicts that roughly 6 per cent of the total 
gas mass is traced by single-component BLAs, which represents 
about 10 per cent of the total WHIM mass in this model; if the 
restriction that the absorbers be single-component is dropped, then 
around 25 per cent of the total gas mass (40 per cent of the WHIM 
mass) can be recovered from the detected broad Hi absorption 
(Sec.|4~5). 

• Although some of the gas mass with temperatures 
T > 5 x 10 4 K and densities A < 10 2 (i.e. the diffuse warm-hot 
phase) can be traced using BLAs, a significant fraction remains 
undetected as a consequence of a minimum (instrumental) sensitiv- 
ity limit. Detection of the bulk of warm-hot gas requires a sensitivity 
(in terms of the Hi central optical depth) of log To < -2 (Sees. [42] 
1431. 
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APPENDIX A: LINE-FITTING 

We fit our spectra using a significantly modified version of autovp 
l |Dave et al.|1997[ l, assuming each absorption component to be de- 
scribed by a Voigt profile given by the analytic approximation of 
Tepper-Garcia (2006). autovp identifies absorption features using 
the equivalent-width significance criterion of Lanzett a et al.| ( |T987) . 
The spectrum is scanned using a window of width n pixels in search 
of regions with significant absorption. A region is considered sig- 
nificant in absorption if its equivalent width satisfies W > N <t w , 
where N is the significance level, and cr w is the uncertainty in the 
equivalent width, integrated over n pixels, given by 



o- w ~ V^|^lol I I +:)-(.V/A'f' 



(Al) 



Here, Ad is the pixel width, S/N is the adopted signal-to-noise ratio, 
z is the (central) redshift of the absorption feature, Aq is the rest- 
frame wavelength of the transition considered (e.g. H i Lya-), and c 
is the speed of light. We adopt N = 7 and n = 25 (corresponding to 
approx. 88 kms 1 for our chosen resolution), and Av = 3.5 km s" 1 . 
With these values, the significance value translates into a rest-frame 
equivalent width 



W r ~ 500 (S/N)" 1 mA. 



(A2) 



Note that our adopted window width does not affect the parameters 
of the fitted line(s) in any way. 

If fitted by a single component, this implies that our line sample 
is formally complete down to H i column densitie^j 



N H , 



: 9 x 10 13 (S/N)" 1 cm" 2 . 



(A3) 



For Doppler parameters in the range £> H i > 40 km s" 1 , characteristic 
of BLAs, the above corresponds to a sensitivity limit in terms of 
absorption strength 



Nm/cm 2 
£>Hi/kms -1 



< 2.3 x 10 12 (S/Nr 



(A4) 



which is equivalent to a Hi Lya- optical depth at the line central 
t > 1 .74 (S/N)" 1 . Note that the values implied by the above equa- 
tion are below the value commonly adopted for the identification of 
BLA candidates in real QSO spectra (see eq.|2j. 

A spectrum is fitted in two steps. In the first step, an absorp- 
tion component is fitted at the pixel with the minimum flux in each 
detection region, starting with the region with the overall minimum 
flux. The column density, Nm, and the Doppler parameter, bm, of 
the line are both iteratively reduced by a factor 0.99 starting from 
large values (e.g. Nm = 10 20 cm" 2 and bm = 300 kms" 1 ; see below) 
until the flux at that pixel is within 2cr below the actual flux level, 
i.e. in the range [F - 2a, F], where cr is the local noise. Further 
components are added and their parameters correspondingly adjus- 
ted, taking all previous fitted lines into account, until the residual 
flux (i.e. the difference between actual flux and model flux) across 
the detection region is below 2cr. This procedure is repeated for 
all detected regions. In a second step, the line parameters (velocity 
centroid vq, column density Nm, and Doppler parameter bm) of all 



15 The quoted value is valid only for absorption lines on the linear part of 
the curve-of-growth, which is the case for the majority of the components 
identified in our synthetic spectra. Also note that we do detect lines with 
column densities (and rest-frame equivalent widths) smaller than the quoted 
values, since a detected region can be fitted by more than one component. 

16 The relation between line-strength and central optical depth for the Lya 
line is given by eq. (53) in Appendix|B| 



lines are simultaneously adjusted using the Levenberg-Marquardt al- 
gorithm ( Levenberg 1944; Marquardt 1963 ) as implemented in Press 



et al. 



1992} until the reduced x 2 -value, i.e. the x 2 -value divided 



by the degrees of freedom, is below kL, = 1.2. If convergence is 
not achieved, the sightline is discarded. We note that the fraction of 
discarded sightlines is vanishingly small, and it amounts to 1/5000 
for our spectra at z = 0.25 and none for our spectra in the range 
0^z<0.5. 

Since we do not take higher order H i Lyman transitions into ac- 
count, saturated Hi Lya lines deserve special attention. A pixel 
is considered saturated if the corresponding flux is below 2cr. 
For (S/N) > 10, this implies that the flux is of the order of, or 
lower than, 2cr min (where cr min = 10" 4 ), which is equivalent to a 
Hi Lya- central optical depth To * 8.52 or a Hi colu mn density 
log(N H i/cm" 2 ) as 13.1 + log(fo H i/kms"') (see eq. B3 I. Assuming 
a Doppler parameter bm = 30km s" 1 (which approximately cor- 
responds to the median b-value of our line sample, see Sec. |3.3[ >, 
this corresponds to log(N H i/cm" 2 ) as 14.5. In order to prevent our 
algorithm from severely underestimating the true column dens- 
ity of such saturated lines, and at the same time to avoid fitting 
lines with unrealistically large column densities along a given 
sightline, during the second fitting step we limit the column dens- 
ity of an individual absorption lin^jto five times this value, e.g. 
log(N H i/cm" 2 ) max as 15.2 for bm = SOkms" 1 . Note that this value 
is not a strict limit but may still vary (in particular, it can be lar- 
ger) depending on the actual fe-value of the saturated line. As a 
consequence, we highly underestimate (by up to two orders of mag- 
nitude; see Table|2j the actual baryon content in Hi along all fitted 
sightlines, which is dominated by high H i column density gas. Yet, 
the properties of the H i absorbers, in particular the broad H i ab- 
sorbers that are relevant for the present study, are not affected, since 
these are dominated by the low column density population with 
log(N Hl /cm- 2 ) < 14.5. 

During the second fitting step, we impose a minimum line 
width of bm = 10 km s _I , corresponding toI« 6000 K assuming 
pure thermal broadening. Lines narrower than this are discarded 
during the fitting process unless the new^ 2 -value increases above 
1.2 k .. We note that this cut-off does not appreciably affect the res- 
ulting b- value distribution (see Fig. [2). Indeed, observations indicate 
that narrow (i.e. bm < 15kms~') Hi absorbers at low redshift are 
scarce ILehner et al.|2007> . 

Since our synthetic spectra are continuum-normalised by con- 
struction, we do not fit a continuum prior to line identification. We 
limit the line width to a maximum value of bm = 300kms _1 , al- 
though larger values are allowed if doing so reduces the x 1 - value 
below 1.2 kLj. As a consequence, we find a small number of very 
broad (bm > 200km s ), very shallow absorption features, which 
become less numerous with decreasing S/N. Since most of these 
features are real, though scarce, we do not discard them but include 
them in our resulting line sample. 

Finally, any candidate lines with formal relative errors in N H , 
or bm larger than 50 per cent are sequentially discarded unless the 
^ 2 -value increases above 1.2^? .. Note that the final formal errors 
are typically much smaller than this, around 10 per cent for both 
N H i and b m - 



17 This particular value chosen is arbitrary, but has been found to give 
satisfactory results. Note that our algorithm may still underestimate the true 
column density of heavily saturated lines. 
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APPENDIX B: OBSERVABILITY OF H i ABSORBING GAS 

The neutral hydrogen column density Nm is given by 
N H i = N H • f H i , 

where N H and f H i = (n H i/n H ) are the total hydrogen column density 
and neutral hydrogen fraction, respectively. The total hydrogen 
column density can be written using the hydrogen particle density 
n H as 

N H = I rig dl = 5h • L , 



fc = f 

Jo 



where L is the physical, linear extension of the absorbing structure 
along the sightline, and iTh is the average hydrogen particle density. 
In the following we will write n H = Sh, but the reader should keep 
the (slight) difference in mind. 

The width of an H i absorbing line as measured by the Doppler 
parameter b Hi may be modelled as 

b 2 m = bl + b 2 H + b 2 „,. 

The thermal width, i.e. the broadening due to the temperature T of 
the absorbing gas is given by 



[fcr/kms -1 ] = 12.9V7710 4 K, 
and the Hubble broadening by (see e.g. Schaye 2001 1 

b H ~ \h(z)-L, 



where H(z) is the Hubble parameter (at the appropriate epoch), 
expressed as 



H(z) = h(z) ■ 10 2 km s _1 Mpc _ 



with 



h( z ) = ho [a m (i+ z ) 3 + n A ] 



1/2 



(Bl) 



(B2) 



We adopt the cosmological parameters {fl m , Qa, ho] = 
{0.238, 0.762, 0.73) as derived from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 3-year data, and find, 
e.g. h(z = 0.25) = 0.81. 

The remaining term, b,„, in the expression for b ai includes all 
other forms of non-thermal broadening and is less straightforward to 
model. It may include turbulence within the absorbing gas, peculiar 
motions of the absorbing structures, etc. Assuming that these are 
negligible compared to the thermal and Hubble components, the 
line width can be approximated by 

[i Hl /kms 4 ] 2 « 166-[r/10 4 K] + ^10 4 [/!(z)-L/Mpc] 2 . 

Putting all the above equations together and simplifying, we find 
that the H i Lyo- absorption strength of the gas is given by 

/ Nm/cnr 2 \ 6.17 x 10 12 h{zT l ■ [n H /10- 5 cnT 3 ] • [WIO -5 ] 



\Wkms l ) ^6.64 x 10~ 2 ■ [77 10 4 K] ■ [h(z) ■ L/Mpc]- 2 + 1 

The central optical depth of the H i hya transition can be expressed 
in terms of (N H i/£>hi) as 



To 



V5e% /Nh, 
m e c \ b m 



■■ 7.56 x 10" 



N H i/cm 2 
bm/ km s -1 



(B3) 



Using the above equations we obtain an expression for t in terms 
of n H , T, and L: 

4.66 h(zT l ■ [n H /10- 5 cm- 3 ] • [WIO" 5 ] 



To 



V6.64 x 10- 2 • [77 10 4 K] • [h(z) ■ L/Mpc]- 2 + 1 



(B4) 
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Figure Bl. H i hya central optical depth, To, as given by eq. (B6J, required 
to detect gas with a given temperature for a range of densities at z = 0.25. 
The value along each curve indicates the corresponding logarithmic hydro- 
gen particle density, log (nn/cm -3 ^. The number in parentheses indicates 
the corresponding logarithmic overdensity. For reference, we include the 
alternative y-axis which shows the corresponding value of the H i Lyo- ab- 
sorption strength at the line centre, log (Nm/^Hi) in units of cirT 2 / km s~' 
(see equation |B3fr . The dashed horizontal line indicates the adopted BLA 
detection threshold given by eq. |2j. 



Note that f Hl = fm(nH, T, z), where the z-dependence comes about 
through the redshift dependence of the ionisation background in- 
cluded in the calculation of f Hl . 

If we assume that the absorbers have linear sizes of the order 
of the local Jeans length (Schaye 2001 ) 

-5 -3 n -l/2 [r/10 4 K] l/2 . [/ g /0.168] 1/2 , 

(B5) 



Lj = 0.169 Mpc [n H /10~ 3 cnT 3 r 1/2 
we get 



4.66 h(z)- [ ■ [n H /lQ- 5 cm- 3 ] ■ [f H ,/l(r 
V2.32 • [n H /10- 5 cm- 3 ] • h{zY 2 + 1 



(B6) 



where we have assumed the fraction of mass in gas to be close 
to its universal valu^j/j, = Q. b /Q.,„ = 0.168. Note that the above 
equation does no longer depend explicitly on the temperature, but it 
does depend implicitly on it through the dependence on fm. 



n H = — *h (1 + z) 3 A - 1.88 x 10-W 3 (-§-) (1 + z) 3 A 
m H V 0.752/ 



Using 

<Pb> 
m H 

the above equations can all be expressed in terms of the overdensity 
A as well. 

Fig. |Bl| shows the Hi Lya central optical depth, To, as a func- 
tion of gas temperature for a range of densities typical of inter- 
galactic gas, as given by eq. l |B6[ >. The alternative y-axis shows 
the corresponding values for log (Nhi/^hi) in units of cirT 2 / km s~' . 
The value along each curve indicates the assumed logarithmic hydro- 
gen particle density, log (n H /cm 4 ), and the value in parentheses in- 
dicates the corresponding logarithmic overdensity, log A, at z = 0.25. 
The horizontal dashed line indicates our adopted sensitivity limit 

18 The total mass density parameter is Cl,„ = CI/, + Cl c . The most recent 
measurements of the baryonic and dark matter density parameters yield, 
respectively, il/, 

[20TT) . 



: 0.0449 ± 0.0028 and Sl c = 0.222 ± 0.026 ( Jarosik et al. 
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Figure CI. CDDF for different model runs in a box L = 100/i~' Mpc per side 
at z = 0.25 using 5000 random sightlines. Dashed lines show the distribution 
of column densities for individual absorption lines obtained from fitting 
the synthetic spectra adopting S/N=50. Solid lines show the distribution 
of integrated Nh i along each sightline. Note that the adopted S/N value is 
higher than the S/N of the data, and the latter are only included for reference. 
The distribution of column densities for individual absorption components 
extends to much lower values since a spectrum (corresponding to a single 
physical sightline) is generally fitted with more than one component. 

as given by eq. l|2j. This figure demonstrates that the detectability 
of H i absorbing gas at a given density drops sharply with temperat- 
ure. For example, the detection of gas with log(nH/cm ) = -4.75 
(which corresponds to A ~ 50 at z = 0.25) and log (77 K) ~ 6 re- 
quires a minimum sensitivity which is two orders of magnitude 
higher than the one required to detect gas with the same density 
and log (77 K) ~ 5. The reason behind the strong dependence of 
log (Nhi/^hi) on temperature is that, at a given density, r is com- 
pletely dominated by the neutral fraction, fm. 



APPENDIX C: CONVERGENCE WITH RESPECT TO 
THE PHYSICAL MODEL 

In this section, we demonstrate that the predicted H i observables 
are robust with respect to variations of the adopted physical model. 

Fig. |Cl| shows the Hi column-density distribution (CDDF) for 
different models run in a box L = 100/T 1 Mpc per side at z = 0.25, 
using 5000 random sightlines. Dashed lines show the distribution 
of column densities for individual absorption components obtained 
from fitting the synthetic spectra adopting S/N=50, as described 
in Sec. [3] solid lines show the distribution of the integrated H i 
column density along each sightline. The latter are important in 
order to remove the uncertainty introduced in the CDDF by our 
fitting algorithm. We include various data sets in this figure for 
reference, but note that our adopted S/N value is higher than the 
S/N of the data. The distribution of column densities for individual 
absorption components extends over a range of much lower values 
since a spectrum (corresponding to a single physical sightline) is 
generally fitted with more than one component. Also, as mentioned 
in Sec. [3] our version of autovp tends to fit (highly) saturated lines 
with more than one component, thus generally yielding lower N H , 
values for each individual component with respect to the Nm value 
obtained by integrating over the corresponding absorption feature. 

Fig. |C2| shows the distribution of Doppler parameters of 
individual components identified in 5000 spectra at z = 0.25 



Figure C2. Distribution of Doppler parameters of individual components 
identified in 5000 spectra at z = 0.25 with S/N=50 adopting different phys- 
ical models in a L = 100ft - ' Mpc box. The dashed vertical line indicates the 
minimum Doppler parameter during the fitting process (see Sec. [3). Note 
that feedback (both by SNe and AGN) leads to a slightly larger number of 
broad lines, although the effect is small. 

with S/N=50 adopting different physical models run in a box 
L = 100/T 1 Mpc per side. Note that feedback (both by SNe and 
AGN) leads to a slightly larger number of broad lines, although the 
effect is not large. 

Finally, Fig. |C3| shows the corresponding fc H i _ Nm distribution 
of individual components. The bm - Nm distribution for each 
model has been obtained as described in Sec. |3.4| For each model, 
the corresponding median values at each bin are connected by 
continuous lines. For simplicity, only the result for model AGN 
displays x- and y-error bars indicating, respectively, the bin size 
and 25-/75-percentiles; but note that the bin-size is identical and 
the scatter similar for all other models. Although the models differ 
slightly from each other, the various Z? H i _ N Ht distributions are 
fully consistent with each other. 

Clearly, feedback by supernovae and AGN does not signific- 
antly affect the observed properties (H i column density distribution, 
Doppler parameter distribution, bm - H i correlation) of the gaseous 
structures giving rise to Hi absorption. This result is consistent 
with previous results by Th euns et al.| j2002) who, however, did not 
consider AGN feedback or metal-line cooling. 



APPENDIX D: NUMERICAL CONVERGENCE 

In this section, we address the convergence of our results with re- 
spect to varying the box size and the mass and spatial resolution. To 
this end, we compare the H i CDDF and Doppler parameter distribu- 
tion of individual components identified in 5000 spectra obtained 
from simulation runs with different box sizes and resolutions, which 
all adopt the model REF at z = 0.25. Note that the use of this par- 
ticular model does not affect our results, since we have shown in 
Appendix [C] that both the Hi CDDF and the distribution of Doppler 
parameters are insensitive to the adopted model. In the following, 
the simulation runs we use are denoted by LxxxNyyy, where xxx 
corresponds to the linear size of the cubic box in h Mpc, and yyy 
to the number of (dark matter, baryonic) particles per side. 

To investigate the convergence with the box size, we use 
the simulation runs L025/V128, L050/V256, and L100/V512, which 



© — RAS, MNRAS 000,[TpO| 



24 Tepper-Garcia et al. 



10' 



10' 



NOSNJMOZCOOL 
NOZCOOL 
REF 
AGN 




10' 



10' 



10' 



10' 



N H | [cm-'] 



Figure C3. b\\ i - Nh i distribution of individual components identified in 
5000 spectra at z = 0.25 with S/N=50 adopting different physical models 
in a box L = 100/i -1 Mpc per side. For clarity, errors bars (corresponding to 
the bin size along the ;t-axis and to the 25-/75-percentiles along the y-axis) 
are shown only for the AGN result. We note that the scatter for the other 
models is similar. 



all have the same mass- and spatial resolution. The convergence 
with resolution is investigated using simulations run in a box 
L = 50/7~' Mpc per side, and varying the (dark matter, baryonic) 
particle number; more specifically, we use the runs L05(W128, 
L05(W256, and L05(W512, whose mass (spatial) resolution varies 
in factors of 8 (2). The choice of this particular box size is arbitrary 
but justified, since our results are converged with respect to the box 
size, as we will show next. 

Fig. |Dl| shows the Hi CDDF and the distribution of Doppler 
parameters using 5000 sightlines obtained from simulation runs 
with different box size but fixed mass- and spatial resolution. The 
top panel shows the distribution obtained from the Hi column 
density integrated along each individual sightline; the middle panel 
shows the distribution obtained from the H i column density of each 
component identified in the corresponding spectra with S/N=50. 
The Hi CDDF of integrated column densities is fully converged 
for column densities N Hr < 10 17 cirr 2 . At higher column densities, 
we do not expect a good convergence, since these column densities 
correspond to the optically thick regime, while our calculations 
assume optically thin gas. The Hi CDDF of individual components 
is fully converged at column densities Nm < 10 15 ctrr 2 , which is the 
relevant column density range for this study. The difference between 
the various Hi CDDFs of individual components in the range 
Nhi > 10 l5 cnT 2 is due to the inability of our fitting algorithm to 
accurately determine the column density of saturated H i absorption 
features. Finally, the bottom panel demonstrates that the Doppler 
parameter distribution in the range 10 km s < £>hi < 100 km s - ' is 
fully converged with respect to the box size. 



Fig. |D2| shows the H i CDDF and the distribution of Doppler 
parameters using 5000 sightlines obtained from simulations run in a 
box L = 50/i~' Mpc per side with different mass and spatial resolu- 
tion. Both the distribution of H i column densities integrated along 
each sightline (top panel) and the distribution of H i column densit- 
ies of individual components identified in the corresponding spectra 
(middle panel) appear to be fully converged with respect to resolu- 
tion at column densities N H i < 10 15 ctrr 2 . At higher column densit- 
ies, neither distribution is fully converged, although the difference 
between the two highest resolution runs, L050/V256 and L05O/V512, 
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Figure Dl. Numerical convergence with respect to the box size at a 
fixed mass and spatial resolution using 5000 sightlines and correspond- 
ing spectra obtained from a simulation which adopts the model REF in a 
L = 100/7" 1 Mpc box per side at z = 0.25. Top: H i CDDF obtained from the 
Hi column density integrated along each sightline. Middle: Hi CDDF of 
the H i absorbing components identified in the corresponding spectra with 
S/N=50. Note the difference in the plotted range between the top and middle 
panels. Bottom: Distribution of Doppler parameters. The dashed vertical line 
indicates the minimum Doppler parameter allowed during the fitting process 
(see AppendixjA). 
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is very small for the integrated column densities. In the case of the 
Hi CDDFs for individual components at Nm > 10 15 cirT 2 , this is 
again in part due to the difficulty in determining the true column 
density of saturated H i absorption features. 

The Doppler parameter distribution shown in the bottom panel 
of Fig. |D2| indicates that the resolution of the L05(W128 run is not 
high enough. The distribution in the L05CW256 run, which has the 
same resolution as our fiducial run, does not show full convergence 
at Doppler parameters in the BLA regime (fo Hl > 40kms~') 
with respect the higher resolution run, 11)50.^512, although the 
difference is small (see also Th euns et al.|1 998 1. 

In summary, our results are robust with respect to varying the 
size of the simulation box, and our adopted resolution is high enough 
to guarantee the convergence of our results in the range of column 
densities relevant for this study. However, the distribution of Doppler 
parameters is slightly sensitive to the adopted resolution in the range 
of interest for BLAs, i.e. for Doppler parameters b ai > 40kms~ 1 . 
Thus, some caution is advised when interpreting results based on or 
making predictions for this observable. 
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Figure D2. Same as Fig. |Dl| for the numerical convergence with respect to 
the mass and spatial resolution at a fixed box size. 
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